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SUMMARY 


The problem of stall flutter is approached in two ways. First, , 
using the M.I.T.-NACA airfoil oscillator, the aerodynamic reactions on 
wings oscillating harmonically in pitch and translation in the stall 
range have been measured, evaluated, and correlated “v^ere possible with 
available published data, with the purpose of providing empirical infor- 
mation where no aerodynamic theoiy exists. The major effects of Reynolds 
number, airfoil shape, and reduced freq.uency on the aerodynamic reactions 
have been reaf finned. No. Instances of negative dan 5 >ing were observed in 
pure translatory motion and the ranges of negative dan^plng occurring in 
pure pitch had the same general trends noted by other experimenters. 

Data on the time-average values in the stall range of both lift and 
moment are presented for the first time. 

Second, the results of numerous experimental observations of stall 
flutter have been reviewed and the various known attempts at its pre- 
diction have been examined, con5>a37ed, and extended. The shaip) drop in 
critical speed and change to a predominantly torsional oscillation 
usually associated with the transition from classical to stall flutter 
is apparently primarily but not entirely caused by the mrked changes 
in moment due to pitch. Fairly good stall- flutter predictions have 
been reported only when adequate en 5 )irical data for this aerodynamic 
reaction happened to be available for the desired airfoil shape, 

Reynolds number range, and reduced- frequency range. A semieB 5 )irical 
method of predicting the variations of moment in pitch with Airfoil 
shape, reduced frequency, initial angle of attack, and an 5 >lltude of 
oscillation has been presented. 


INTRODUCTION 


Stall flutter differs from classical flutter in that the flow over 
the airfoil is stalled during all or part of the cycle of oscillation. 
Althou^ stall flutter has been observed in aircraft propellers for 
some time, new interest in the phenomenon has been generated by the 
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advent of the aircraft gas turbine -vdiere turhine-blade stall flutter 
is a serious problem. Even thou^ it is known that the critical speed 
drops sharply from the value predicted by classical theory as portions 
of the flow begin to stalls no adequate prediction techniques have yet 
been proposed primarily because not enou^ is known about the behavior 
of the aerodyneimic reactions in the osci3JLating partially stalled flow. 

Work on staUL flutter t/as begun at M.I.T. in the simnner of 19^7 
when the evaluation of the low-angle-6f -attack aerodynamic-derivative 
data was conpleted. The M.I.T.-NACA airfoil oscillator was modified 
to operate in the stall range and a large mass of aerodynamic-derivative 
data was accumulated which within its someidiat limited range gives a 
far more conplete plct\ire than has been heretofore available. 

In this report an attenpt has been made to evaluate these data with 
conparison^ where possible, with other available data. The information 
thus obtained iras used in critically examining proposed stall- flutter 
prediction techniques and actual stall- flutter tests. 

This work was done at the Massachusetts Institute of Technology 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 


0) 

<1 

a 

h 

P 


b 


angular frequency of forced motion 


dynamic pressure 



pitching angle of wing, positive in direction of stall 


vertical translation of wing at 37-percent chord, positive 
downward 

angle between front and. rear actuator wheels or, as defined 
in reference 1, angle between pitching and translatory 
motion 


semi chord 


V air-stream velocity 

K reduced- frequency parameter (o5b/V) 

0 angle by which pitching motion leads translational motion 
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Description of Apparatus 

Actuator and lin'kage .- The test-section arrangement is shown in 
figure 1. The mechanical features of the oscillator are shown diagram- ■ 
matically in figures 2(a) and 2(b)^ which also briefly outline the func- 
tions of the various parts. The numbers in parentheses in the descrip- 
tion below correspond to numbers of figures 2(a) and 2(b) . 

Power is supplied to the line shaft (l) by a 15-horsepower electric 
motor with a Ward-Leonard speed control. V-belts drive sheaves ( 2 ), 
which are geared to drive the crank wheels (3). These crank -vdieels all 
have variable-anplitude cranks and connecting rods which drive the 
vertical cross heads (i^-) at oscillator frequency. The phase P between 
the front and rear pairs of crank -vdieels may be varied in increments of • 
3 . 75 ° with the use of the angular scale (15)- The position indicator (6) 
gives the position of the rear vertical members and the position of the 
other members may be found if the phase angle p between the front and 
rear crank--ydieel pairs is known. 

From the vertioal cross heads the motion is carried to the cam ( 16 ) 
by steel bands (7). These bands are attached to the vertical cross ' 
heads by tumbuckles and to the cam by forks operating in ball bearings. 

A similar set of bands connects the cam to the springs ( 9 ) in the over- 
head stincture, and the req.uired initial band tension is attained by 
use of electric motors and gear trains (lO). The motion of the cam is 
transmitted to the wing and accelerometers (I 7 ) throu^ a four-bar 
linkage. The vertical bar ( 18 ) at the cam center is free to move verti- 
cally so that the cam may tianslate vertically, pitch, or move in any ‘ 
combination of the two motions. 

The three airfoils used in these tests were constructed of sycamore 
wood and 0.007- inch magnesium-alloy sheet. The uncovered frame is shown 
in figure 3* Airfoil section and ordinates are given in figure ^4-. 

Instrumentation . - The instrumentation used in these tests is that 
of flgui*e 5* The amplifying equipment, the bridge balances, and the 
recording oscillograph are products of the Consolidated Engineering 
Corporation. 

( 1 ) The strain gages which measure lift, moment, drag, position, 
and translational and pitching accelerations are supplied with 1000-cycle 
carrier voltage by an oscillator nanufactured by the Consolidated 
Engineering Corporation. Four strain-gage pairs are used to measure 
lift, two for drag, two for moment, and two each for pitching and trans- 
lational acceleration. 

(2) The equalizer panel consists of 15-ohm potentiometers in series 
with each strain gage (200 ohms) so that the sensitivity of each strain 
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gage may "be varied Indlvldiially; thus the strain- gage outputs may he 
made equal under a given load. 

(3) The offset circuit utilizes resistances to unbalance the 
bridges by a known amount. Its purpose is to balance out the poortlon 
of the strain-gage signal due to static lift, drag^ or moment, at hi^ 
angles of attack. Thus the oscillating portions may be rendered on the 
sensitive paper of the oscillograph at a reasonable aji 5 )litude, keeping 
the anqplifler output within the linear range. 

(4) The bridge balance units consist of both a phase and resistance 
balance and , are used to aidjust the strain-gage output to zero under any 
tare loading. 

(5) The position indicator serves to indicate reference points in 
the wing’s cycle of motion. 

(6) The attenuator boxes are used only in the lift and moment 
Inertia circuits for balancing the inertia signal from the acceleromr- 
eters against that of the wing. A better understanding of their use 
n.n(^ operation may be gained by reading the section on calib 2 rations. 

(7) The anpllflers, as mi^t be supposed, anplify and rectify the 
modulated 1000-cycle signals from the strain gages, for use in the 
oscillograph. The injected carrier shown feeding into the an 5 )lifiers 
in figure 6 is used to provide a reference signal an 5 )lltude level so 
that the anpllfier output shows the positive or negative sense of the 
signal. 

(8) The add circuit box serves the purpose of subtracting the 
inertia signal (from accelerometers) from the total (aerodynamic plus 
inertia) signal in lift and/or moment. The box is equipped with a 
switch so that the operator can control whether the signals shall pass 
through the add box and appear separately on the oscillograph record, 
or -vihether the subtraction described above shall take place and only 
the aerodynamic portion of the signal be recorded on the sensitive 
paper. 

(9) The recording oscillograph utilizes hl^-natural- frequency 
galvanometers which deflect a beam of li^t onto sensitive paper for 
recording puiposes. There are l4 channels available for use, 8 Of 
which were used in these tests. The quantities measured together with 
their trace numbers are listed below: 


1 

Lift 

9 

Lift inertia 

3 

Drag 

10 

Moment inertia 

5 

Position 

u 

Lift inertia 

7 

Moment 

12 

Moment inertia 
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D 

L 

M 

0 

C 

Tr 

Kr 

acm 

X 

F 

a 

C(K) 

a 

c 

Ah 

Ax 

t 

6a 

>^3 

Xc 


drag 

lift force on airfoil 
moment on' airfoil 

phase angle hetwsen motion and reaction for sinusoidal 
reactions or lag angle as defined in reference 1 

coefficient 

characteristic time, defined in section, "Stall-Flutter Analysis" 

reduced frequency Invol-vlng characteristic time Tj. 

varlatlm of from static curre ‘ <«steaay-state) 

dynamic slope of moment against angle- of-at-tack curve 

Theodor sen *8 function as defined in reference 2 (F(K)) 

Theodor sen’s function as defined in reference 2 (G(K)) 

Theodorsen’s function as defined in reference 2 (F(K) + iG(K)) 

ratio of distance of elastic axis behind midchord point to 
semi chord 

chord of wing 

increment of translatory motion 

incremental angle of attack or Interval between mean angle 
of oscillation and steady-state stalling angle 

time 

interval between dynamic and steady- state s'tall 
^“•stall " “■steady- state stall) 

torsional daji 5 )ing parameter as defined in reference 3 

frequency parameter used in British nomenclature = 2 k) 

dimensionless center- of-gravlty dis'tance as defined in 
reference 2 
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moment of inertia per unit length of wing about a 

static unbalance parameter for wi n g as defined, in reference 2 

dimensionless radius of gyration of wing as defined in 
reference 2 

torsional stiffness parameter as defined in reference 2 
bending stiffness parameter as defined in reference 2 
natural torsional frequency 
natural bending frequency 

mass ratio for wing as defined in reference 2 
work 
Subscripts : 
c critical 

t torsion 

L lift 

M moment 

T due to translational motion 

P due to pitching motion 

i initial or mean geometric position of airfoil 

o an^l.1 tude of fundamental harmonic of motions or reactions 

(e.g. , a = ttoei^) 

M.I.T. EXPKEIIMEN!EAL PEOGEAM 


Sa 

% 

Ca 

Ch 

obc 

CDb 

m/jtpb^ 

W 


The apparatus used in the over -all research (reference 4) was 
modified to make the apparatus -useful for studies at hi^ angles of 
attack. A review of the varlo-us elements of the mechanism both 
mechanical and electrical^ and a con^ilete outline of calibration, 
testing, and analysis procedures is contained in -fche following pages 
of this section. 
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It Td.ll "be noted that lift inertia and moment inertia are each given 
tTTO oscillograph channels . The reason behind this apparent duplication 
is that one lift inertia and one moment inertia trace Tdll not appear 
on the record after subtraction from the appropriate total aerodynamic 
reaction signals. 


Testing Procedure and Analysis of Data, 

Linkage setting and adjustment .- Immediately after assembly of the 
oscillator mechanism in the Tdnd tunnel, careful adjustments of the 
linkage -vrere made. These adjustments consisted of the following: 

( 1 ) Adjustment of an^litude of each crank to prevent any relative 

motion betrreen the left and rl^t wing supports 

( 2 ) Adjustment of tumbuckles to level all elements in the zero- 

aji 5 )lltude position 

( 3 ) Check and adjustment of the clearances bet-ween the Tdng and 

end plates 

After the linkage T/as adjusted, it "was set to give a pure pitching 
motion. This -was acconplished by setting the amplitudes of the cranks 
equal, and with l 80 ° phase angle between the front and rear cranks. 

Phase angles betrreen front and rear cranks were obtained by pinning the 
rear cranks in position, unmeshing the driving gears and setting the 
phase angle Trith the aid of the scale adjacent to the front crank Tdieel, 
and then remeshlng the driTrlng gears. 

After the linkage Treis adjusted to give the desired motion, it -was 
necessary to be able to set accurately the Initial angle. A protractor 
incorporating a very sensitive level was used to measure the angle of 
the Td.ng-support strain-gage beams with the horizontal. The angular 
measurements made "with this device were quite accurate, being Td-thin 
±0.1° of the desired value. In order to be certain that all four "wlng- 
support points were in a plane, a leveling plate was used. After the 
initial angle Tras set on the left side, the ri^t rear element -was 
lowered and the leveling plate was placed over the remaining three 
points. An adjusting screw was utilized to bring the rear element up 
into contact Tdth the plate, the contact being indicated by a deflection 
of the needle on the dial of the amplifier attached to the moment bridge. 
This method permitted very accurate setting of the four supporting 
points so that their' deviation from a plane -was not greater than about 
0.002 inch. 

i-Then the initial angle was set, the strain-gage beams 'were loosened 
and removed, and the wing -was placed between the end plates; then the 
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strain- gage 'beams -were inserted in their slots and tigjhtened^ the 
remova'ble fairing plates replaced, and the mechanical appamtus was 
ready for tests. 

Calihration . - Before any tests were actually run it was necessary 
to cali'brate the equipment. The accelerometer circuits were provided 
with attenuators in the form of voltage dividers so that the accelerom- 
eter output could "be varied from the full output to zero in 0 . 01 -percent 
steps. Since the accelerometers were to he used, to cancel the inertia 
signals of the wing, a device putting out the same Inertia signal as 
the wing was indicated. The desired Inertia signal was provided hy the 
use of a "dummy" wing. The dummy was placed in the linkage in the same 
nmtnner as the win^ and oscillated in a desired motion. The signals 
from the accelerometer and the dummy wing were then added and fed into 
a cathode-ray oscilloscope. By varying the output of the accelerometer 
with the voltage divider described above it was possible to cancel out 
almost exactly the dummy-wing Inertia signal. The "best apparent setting 
of the voltage divider was chosen and the added signals were fed into 
the recording oscillograph. Records were taken at the best apparent 
voltage-divider setting and at several values on either side. After 
these oscillograph recordq were developed it was possible to find one 
in which the resultant signal was zero. The voltage-divider setting 
corresponding to this near perfect cancellation was then used throu^- 
out the entire series of tests. 

If an indication of the average force or moment on a unit span of 
the wing was to be measured, it was necessary to have the sensitivities 
of all strain gages measuring that quantity equal. This was accomplished 
by applying an equal load to each supporting beam in turn, .and adjusting 
the appropriate potentiometers in the equalizer box until the output for 
any .one of the positions of loading was the same as for the others. 

This procedure was carried out for each of the lift, moment, and drag 
measuring elements. 

Having equalized the sensitivities of the force measuring elements, 
it was possible to calibrate. For lift, this was accong)liBhed by placing 
a wooden platform over the dummy wing and loading it with wel^ts of 
known value. Moment was calibrated by hanging known wel^ts on a bar 
at the center line of the rear elements, while drag calibration required 
a special rig of wire and pulleys so that loads could be applied in a 
horizontal plane. 

In a]JL these calibrations tare loads were taken into account and 
zero checks were made after every series of loads. A complete calibra- 
tion of force measuring elements was made after every series of tests, 
and at all values of attenuation used during the tests. 
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A previous tunnel- throat survey indicates negligible .variations in 
velocities across the test section; therefore no time Tvas spent in making 
a similar survey^ since conditions in general vere almost Identical. The 
tunnel calibrations were limited to a measurement of the velocity at the 
wing position in coirparison with the velocity at the tunnel .pi tot- tube 
position. 

¥ind-on tests .- The first wind-on tests were conducted to obtain the 
static airfoil characteristics of the three wing sections. These tests 
were run at 95 inlles per hour, and consisted of a measurement of the lift 
and moment on the airfoils at intervals of angle of attack of approxi- 
mately 3° from the angle of zero lift up throu^ the stall angle. The 
results of these tests ara plotted in figures 8, and 9* 

It was found that three men were req.ulred to run the systematic tests 
One to operate the wind tunnel and oscillator, one to operate the instru- 
mentation, and one to take data and check the steps of the instrumenta- 
tion operator. Tests were- carried out throu^ the conplete oscillator 
frequency range on each of the three airfoils at initial angles of 0°, 

6°, 10°, 12°, ll)-°, 18°, and 22° in pitch and translation. 

The test procedure was about the same for all of the tests. With 
the wings Installed between the end plates, and at the proper mean angle 
of attack, all necessary electrical adjustments were made and the tunnel 
was brou^t up to speed. 

With the manometer set at if. 00 inches of alcohol the static lift 
force and moment signals were offset to as near zero-airspeed position 
as possible, and the oscillator was started and run up to the desired 
test frequency. 

With the oscillator at test frequency the Inertia and aerodynamic 
signals 'vrere adjusted separately to be just within the linear range of 
the anpliflers. The signals were then added and a record was taken. 

The above procedure was repeated for each frequency, and then the 
tunnel power was shut off. 

Men the tunnel airspeed had reached zero, the offsets were all 
adjusted to the no-load position, and zero recoi^is were taken at all 
values of attenuation used in the tests. 

Readings of tenperature were taken for each series of tests so that 
necessary corrections in velocities could be made. 
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IThe procedure followed in each test is outlined again in the 
following check-off list: 

(1) Check carrier voltage (lO volt) 

(2) Balance bridges (injected carrier "off") 

(3) Set injected carrier switch to "strain gages" 

(4) Switch add box on "sep" 

(5) Locate galvanometer zeros for lift, moment, and drag channels 

(6) Put oscillator in neutral (pitch only) 

(7) Start tunnel and bring up to speed 

(8) Offset lift, moment, and drag to galvanometer zero positions 

(9) Start and bring oscillator up to test frequency 

(10) Set attenuator for maximum allowable an^litude (all traces) 

(11) Set all meter switches on "galvanometer" 

(12) Switch add box on "add" , 

(13) Take a record 

(14) Repeat steps (8) throu^ (13) for each frequency 

Systematic tests .- The two types of motion considered are: (l) A 

pure pitching motion about the 37-pe^6^'t chord with an amplitude of 
6.08° and (2) pure translation with an anplltude of O.9 inch. 

The first series of tests in both pitch and translation were -carried 
out at about 95 miles per hoxir and Included tests at singles of 0 , 6 , 

10°, ll4-°, and 18° with frequencies of 4, 8, 12, and 16 cycles per second 
for each airfoil at each angle of attack. A supplementary series of 
tests was devoted entirely to pure pitching motion with the threefold 
purpose of extending the range of K, showing trends at critical angles 
of attack not included in previous tests, and showing that small changes 
in Reynolds number would not affect the data to a large extent. 

The sipplementary tests were ceirried out at about 95 miles per hour 
except where noted and Included (l) tests of each airfoil at each of the 
angles of attack listed above at frequencies of 2 and I8 cycles per second; 
(2) tests on each of the airfoils at mean angles of attack of 12°, l6°, 
and 22°, and (3) tests on the Intermediate airfoil at a mean angle of 
attack of l4° at frequencies of 2, h, 8, 12, I6, and I8 cycles per second, 
with velocities of 65, 80, 95, and 110 miles per hour. 

Record analysis .- The analysis of records described below consists 
of a measurement of peak amplitudes and apparent phase a ng les and a 
harmonic analysis to determine net work per cycle and fundamental ampli- 
tude and phase of moment in pure pitch. The oscillograph traces utilized 
in this analysis are lift, moment, and position. No data on drag are 
presented in this report since a preliminary analysis of drag data indi- 
cated a coupling of lift and drag in such a manner that at times a con- 
siderable error in drag was Introduced. It appears, however, that drag 
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forces are of little significance in the hl^-angle-of-attack flutter 
prohlem, so that this omission is relatively unimportant.- 

The lift and moment Inertia traces were not used directly in record 
analysis. The presence of these traces on the record served as a con- 
tinuous check on general accelerometer performance and also provided a 
basis for spot checks 'vAilch con5)ared the measured accelerations with 
calculated values. 

The lift and moment traces are representations of the oscillating 
aerodynamic lift and moment, so that measuring the an^lltude of the 
lift or moment trace and multiplying it by an over-all calibration 
factor give values of oscillatory lift or moment in pounds or inch- 
pounds, respectively. The coefficients of lift and moment (which are 
plotted in figs. I6 throu^ I9) were calculated by use of the formulas: 


Cl 


L 

4qb(span) 


i^qb^( span) 


The symbols Cj^ and Cj^ were used since they represent a magnitude 

without reference to wave form, whereas \/r^ + which was used in 

reference 4 indicates the an^jlitude of a sinusoid. 

The phase angle ^ is defined as the angle by which the force or 
moment leads the motion. This definition, however, is applicable only 
to a sinusoid. At large angles of attack, both the lift and moment 
traces tend to deviate from the pure sinusoids predicted by flutter 
theory. In cases -rfiere the wave form of the trace is not sinusoidal 
at the peak, no phase angles have been measured. In all cases, at angles 
of attack above 6°, caution must be exercised in the Interpretation of 
measured phase-angle data. The plotted values of 0 are averages of 
data taken from four cycles of each oscillation. 

The positive directions of reactions, the peaks ^associated with 
maximum reactions, the location of peak angles of attack, and the 
distances measured to obtain phase angles are described in the diagram 
of figure 10. 
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Ttie work per cycle for a pure sinusoidal motion and reaction may 
"be determined as follows: 


-f 


Wpitch = da = nOoMosin ^ 


^translation ='0L ^ = «tioLo sin 


i 


Since in general M and L were not pure sinusoids it was desir- 
able to resort to a grapbical integration to obtain work per cycle. One 
cycle of eacb record was chosen and Integrated, and occasional spot 
checks were nade to Insure that a typical cycle was chosen. The incre- 
ments ^ and /^h were related to equal time increments in the nanner 
illustrated by the table of figure 11, -which also shows "the graphical 
layout for the integration. 


The aod-s directions of figure 10 determine the sign of the work so 
that a positi-ve sign indicates work done on the -wing by "the air streeim, 
and a negative sign Indicates work done on the air stream by the -wing. 

Since work per cycle is proportional to -fche magnitude of the 
imaginary conponent of the fundamental harmonic of moment, a similar 
integration from a plot , of moment against a motion 90 ° out of phase 
with the actual motion gives an area -vdiich is proportional to the real 
conponent. The ratio of work per 'cycle to this second area in work 
units is the -tangent of the phase angle of the fundamen-tal. From the 
equation above it can be seen -that: 


Mo = - _ 

Jtao sin ^ 

-where is the an^jlitude of the fundamen-fcal harmonic of the moment. 

Thus an^jlitude and phase angle of the fundamental harmonic of moment 
can be calculated. 

The computation of the time-average values of lift and moment 
involves both measurements from the oscillograph records and the electri- 
cal offsets. First the dis-teince from -the mean line of the lift (or 
moment) trace to the zero-airspeed trace was measured. This dis-fcance 
was converted to a force (or moment) by a calibration factor for -the 
particular attenuation used. The recorded value of electrical-offset in 
terms of force (or moment) -was added to the -value obtained by measurement 
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of distance on the oscillograph record to give the time-average value 
of the reaction. 


AMLYSIS AMD CORRELATION OP EXPERIMENTAL DATA 


The experimental data discussed in the following pages of this 
section are the results of two different types of tests: (l) A series 

of flutter tests on wings of various properties held at hi^ angles of 
attack in an air stream and (2) tests involving the use of mecheinical 
apparatus irfiich is designed to inpart a particular motion to an airfoil^ 
in a wind stream. 

The tests in the category of item (l) above have mostly been con- 
ducted in Europe, thou^ two sets of American dhta are referenced later 
in this section. The data from such tests serve to describe the 
mechanisms producing flutter at hi^ angles of attack in a very general 
way, so that forces and moments producing the oscillations cannot be 
directly related to the motions they produce. In order to relate 
reactions and motions the apparatus of item (2) above has been developed 
together with approprtate electrical measuring equipment. 

Thus the aerodynamic reactions associated with a given motion may 
be accurately determined. The apparatus used in the M.I.T. tests is 
described under Description of Apparatus' of this report and the equip- 
ment used in the British tests is descrtbed in reference 

Lift and moment in both pitching and translational motions are 
measured to obtain the M. I.T. data presented below. The British data 
Involve only a measurement of moment in pitch or the out-of-phase com- 
ponent of moment in pitch. 

The M.I.T. data have been taken with the value of Reynolds number 
near 1 x 10^ while other investigators have used Reynolds numbers between 
1.32 X 1q 5 and x 10^ with the majortty of the data from tests in a 
range below Reynolds numbers equal to 4 x 105. 

Only very hi^ Reynolds number data are applicable to aircraft wings, 
but the data taken with Reynolds number 1 x 10^ and below are in the 
proper range for aircraft propeller and turbine blades. 


M.I.T. Data 

The oscillograph records from 150 tests have been- analyzed eind are 
presented in graphical form in figures 12 throu^ 48 , and in tabular form 
in tables I to IX. 
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The majority of the graphs use K = — as the ahscissa. In most 

of the tests, airspeed has heen held constant so that K is proportional 
to frequency. In cases -where velocity has heen -raided, the velocities 
used are designated on the graph. 


All da-ta presented are -within 5 percent of the true value. The 
maxi mum error in check calihratlons is 2 percent and the instrumentation 
introduces an error no greater than 3 percent. These variations -were 
allowed because of the nature of the da-ta and the extreme amount of time 
and expense required to main-tain greater accuracy with the equipment used 
in -taking data. Some da-ta in the -yicinl ty of the stall are not repro- 
ducible to the accuracy described above because the aerodynamic situa- 
tion here is some-what akin to a state of neutral equilibrium and the 
condition of -the flow- depends on tunnel flow inconsistencies as well as 
on the motion of the airfoil. These data, however, represent a very 
small part of "the total, and data from angles near the stall my be 
accepted as represen-fcative of the situation existing unless it is other-' 
-wise noted on the graph or -babulation. 

Some data have been obtained by a direct measurement from the oscil- 
lograph records of an 5 )litudes and phase angles. Since the reaction 
traces are not sinusoids in mny cases the measured values lose some 
significance. However, the basic trends of the data are well-presented 
in this manner, and the method permits -the analyst to cover a great deal 
of territory quickly. 

Dimensionless coefficients of oscillatory force and moment and the 
phase relation of each of these -with the motion have been measured 
directly from -the oscillograph records and are plotted in figures 12 
throu^ 19* It -will be noted that for low angles of attack the escperi- 
men-tal points con^iare quite -well -wi-th the theory. At hl^ angles of 
attack, however, the traces of force and moment are no longer sinusoids, 
so that direct measurement on the oscillograph records gives ralues 
which represent -the peak ralue of an arbitrary -wave form and the time 
interval between its peak and the peak of the motion. The hl^-angle- 
of-at-tack values are plotted as phase angle and aii 5 )lltude because the 
phase angles particularly have been found to be very close to the values 
obtained from harmonic analysis even -where the trace of "the reaction is 
far from sinusoidal. Further, both phase angles and aii 5 >lltudes have 
been found to have -the same general variation -with K as -fche values 
calculated by the harmonic-analysis method -which involves considerably 
more time and expense. 

Some points in -tlie figures mentioned above have no curves faired 
throu^ them. This nay indicate that the points are from low-angle- of- 
attack data -where the reactions are essentially sinusoidal and hence 
should conpare -well -wl-th -tiieory. Also in some instances at hi^ angles 
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of attack, experimental points have fallen very near theoretical curves 
and have heen left without curves faired through them to avoid confusion. 

The plots of work per cycle (figs. 20 throu^ 26) an d the plots of 
j)hase angles and an^jlitudes from harmonic analysis (figs. 27 throu^ 32) 
are the results of the grapMcal integration described in the section on 
record analysis. 

In the plots of work per cycle, the positive-work srange indicates 
that work is done on the wing hy the air stream, while the negative-work 
range indicates that work is done on the air stream hy the wing. It can 
he seen that positive work corresponds to negative damping so that in a 
free oscillation the amplitude would tend to increase if no other daoiplng 
were present. 

In the plots representing work per cycle in pure pitch the curves 
appear to describe a range of K in -vdiich positive work may he done. 

The graphs of figures 20 through 22 show the lower limit of this positive- 
work range and indicate that an upper limit exists hy their tendency to 
assume a negative slope at hi^er values of K. The data cougpiled hy 
Victoiy in refei^ce 3 and discussed in the following section show a 
continuation of this trend at values of K hi^er than those plotted 
in figures 20 through 22. 

In the p'lots representing work per cycle in pure translation, the 
curves all tend to remain in the negative-work region. At* very low 
values of K (much lower than were obtainable with the equipment avail- 
able in this program) it appears that there may he a possibility of posi- 
tive work. Such a possibility is discussed by Den Hartog in reference 6 
in connection with transmission line galloping. A further discussion of 
this possibility is contained in the section Stall-Flutter Analysis. 

The graphical integration described in the section on record analysis 
has been used to determine the fundamental aB5)litude and phase angle in 
the trace of moment due to pitching. This type of analysis parallels the 
analysis of reference so that it is possible to cou^are the trends of 
the two sets of data (see following section). Further, the trace of 
moment in pitch, Tdille periodic, is farther from a pure sinusoid than any 
of the other reaction traces; and the fact that there is good agreement 
in the general trends of the data from the harmonic analysis and data 
measured directly from oscillograph records serves as an indication that 
the measured data (figs. 12 throu^ 19) nay be accepted as representing 
the trends to be expected from a harmonic analysis of the same basic 
data. A conparlson of the phase angles of figure 12 with those of fig- 
ures .27, 28, and 29 shows quite a remarkable agreement in order of magni- 
tude and slope of the curves. 
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In general the coefficients of moment in pitch, plotted from the 
data of the harmonic analysis ^fig. 3P)» are lower in value than the 
coefficients measured directly from the oscillograph records (fig. l6) . 
This is to he expected since there are obviously hi^er harmonics which 
will contribute to the peaA anplitude. 

Figures 31 and 32 are respectively phase angle and amplitude of 
the fundamental harmonic of moment in pure pitch for various airspeeds. 
These tests were made on the intermediate airfoil to find ' if the test 
data were severely affected by Eeynolds number at hi^ angles of attack. 

A mean angle of attack of lk° was used in all tests. This is the stalling 
angle of the intermediate airfoil at 95 miles per hour. 

Both of the above plates show some effects due to velocity changes. 
However, the trends of all curves are similar, and it is suspected that 
the primry reason for the differences is the variation of stalling angle 
with Reynolds number. 

Figures 33 throu^ 36 are presented to indicate the possible values 
of time-average lifts or moments under oscillatory conditions and to 
conpare these with steady-state values. These coefficients could be 
used to calculate the total mean force or moment on the airfoils under 
oscillatory conditions. 

Figures 33 and 3 !*- are plots of the time-average moment coefficients. 
Figure 33 shows the time-average moment coefficients due to pitching. 

These vary considerably from the steady-state coefficients even at mean 
angles well below the stall. This deviation occurs on all of the curves, 
so that it is too consistent to be attributed to experimental error. It 
seems that this variation must be- due to a change in center of pressure 
caused by s-talling of the flow -where local angles of at-tack induced by 
the pitching motion are hl^. It is of Interest to note that at l8° mean 
angles the curves of fi^re 33 show that moment increases -with K, while 
at 22° the opposite is true. The curves of figure 3^ show only relatively 
small de-vlations from the steady-s-fcate curves, and in general it seems 
that the s-balling angles are almost the same -with translatory oscilla- 
tions as in the steady- s-tate case,. At hl^er values of K, the maximum 
-value of lift coefficient is sli^tly greater than the steady- state maxi- 
mum -value in all cases. 

Figures 35 and 36 are plots of the time-average lift coefficients 
in both pitching and translation. The da-ta for the blunt -wing in fig- 
ure 35 are somewhat inconplete since -the zero records for 6° and 10° 
were found in error for all -values of K except 0 . 06 . However, -with 
■the data from the other -two airfoils and from both motions, these data 
seem to indicate that time-average -values of lift below the static stall 
may be assumed to be equal to the steady-s-tate -values without serious 
error. Near -the steady-state stalling angle -the time-average values of 
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lift may either exceed or he less than steady-state values for a pitching 
motion (fig. 35) ^ lu.! the lift continues near the maxi imam value for angles 
well ’above the steady-state stalling angle. The relation of lift coeffi- 
cient to K in figure 35 at hi^ angles of attack is interesting in com- 
parison with the relation of moment coefficient to K in figure 33- As 
in figure 33^ at 18° mean angle the coefficient increases with K, but 
at 22® mean angle, ^ere the variation of moment coefficient with K is 
reversed, the variation of lift coefficient with K remains the same 
thus indicating that between l8° and 22° mean angles the center of pres- 
sure has shifted. , . 

The curves of figure 36 for the blunt and sharp wings indicate that 
the lift in translation is some'vdiat hi^er than the steady-state lift, 
while the data for the intermediate wing shoW' about the same values for 
each. The lift coefficient increases with K in a3JL cases, and the 
maximum values of lift coefficient are hi^er than those in the steady- 
state case, althou^ not necessarily occurring at the steady-state 
stalling angle. 

All of the figures presented above give the results of various 
operations on the basic experimental data, and while they show correctly 
the trends and magnitudes of the quantities involved it is believed that 
the reader may gain a better physical feeling for the problem and a 
better understanding of 'vdiat actually occurs during a cycle of oscilla- 
tion if information such as that in figures 37 throu^ 4l and 46 through 
48 is presented hiia. 

Figures 37 through 4l illustrate the variation of moment in pure 
pitch throughout a cycle of motion. The anplitude of the sinusoidal 
pitching motion is 6.08° in all cases and the pitching axis is at 
37-percent chord. 

Figure 37 illustrates the change with increasing mean angle of 
attack at a constant frequency. It will be noted that at 0° and 6° mean 
angles the moment traces a near elliptical shape quite similar to the 
path predicted in conventional flutter theory. The difference between 
the plotted figures and the elliptical path predicted by theory is due 
in part to the poresence of hash on the oscillogra;^h records and in part 
to the fact that the transcribed points were connected by stral^t lines 
rather theila faired curves. 

At 10° mean angle there is a definite break away from the ellipse 
at the ma-)Hmmn angle of attack, tending to reduce the negative-work area 
enclosed in the loop. At l4 mean angle the top portion of the loop is 
similar in position and curvature to the ellipse predicted by theory, 
but the airfoil stalls sharply near the maximum angle of attack so that 
the loop is enclosed by a clockwise progression of the moment, eind the 
enclosed area is thus a measure of positive work. It is this type of 
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moment vuriation ■whicli causes negative dancing at hi^ mean angles of 
attack. At l8° mean angle the airfoil is always partially stalled as 
is indicated hy the fact that the moment never reaches its steady- state 
value. At 22° mean angle the airfoil is almost entirely stalled through- 
out the cycle and the work area has €igain hecome negative. The varia- 
tions in moment are prohahly due principally to shifts in the center of 
pressure. The negative moment (figs. 37(e) and 37(f)) does in fact indi- 
cate that the center of pressure has moved hehlnd the 37° chord line. 

The small positive-woik area near the m-iniTmim angle of attack indicates 
that the airfoil partially recovers ffcaa the stalled condition, and then 
stalls again almost immediately. 

Figure 38 shows the increase in enclosed area with Increasing 
freq.uency for a pure pitching oscillation about a mean angle of 6°, all 
other conditions being held constant. This effect is predicted by flutter 
theory which is a good approximation to the actual situation at this low 
angle of attack. 

The loops described by the variation of the moment throu^ a cycle 
assume various shapes depending on airfoil shape, mean angle, amplitude 
of motion, and K. Figure 39 illustrates a variation of the frequency 
through the range previously described. Since this oscillation is 
partially in the stall the effect is different from the low- mean-angle 
case described above. It is evident from an examination of these curves 
that the paths are con^osed of a part in which the moment attempts to 
follow the low-angle- of- attack slope of the moment curve, and another 
portion in idilch stall and recovery occur. For low frequencies the 
recovery is made before the ml n .1 mum angle is reached, but at the highest 
frequency con 5 >lete recovery from the stall does not occur until after 
the minimum angle is passed. Figure iK) shows the same general variation 
with frequency as figures 38 and 39/ S’ different mean angle. 

Figure kl ia Included to show types of moment variations not shovm 
in other plates. It can be seen that in general (or at least in the 
case of thesd three airfoils) with increasing mean angle of attack the 
moment through a cycle describes (l) an ellipse as predicted by theory, 

(2) a near ellipse except for a concavity on its upper side, (3) a flgure- 
ei^t-shaped cuive in idilch one loop indicates negative work and the 
other indicates positive work, (4) a loop which indicates all positive 
work, (5) another figure ei^t similar to curve ( 3 )/ and (6) a loop 
which indicates all negative work when the airfoil is at a mean angle 
hi^ enougjh that it never recovers from the stalled condition. 

Several investigators, notably Reid in reference 8, have noted that 
in a pitching oscillation the stall may occur at an angle of attack con- 
siderably above the static stalling single. The amount which the oscil- 
latory stalling angle exceeds the static stalling angle is dependent on 
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a number of factors among -wbich are K, da/dt, airfoil shape, and 
an^jlitude of oscillation. Figures h2, 43, and 44 are graphs showing 
the variation of 6a = - ag^^dy- state stall ^ 

three airfoils of the tests at various mean angles and for an amplitude 

of oscillation of 6.08°. Figure 45 shows plots of: 6a- against a = ^ 

dt 

at stall for the same conditions. On all of these plots the curves are 
faired throu^ the points in order of Increasing frequency and the upper 
points at d = 0 are calculated on the assumption that at a sufficiently 
hi^ frequency the airfoil will not stall before its maximum pitching 
an^jlitude is reached. The data required for all of the 6a plots were 
obtained by measuring the oscillatory stalling angles from plots similar 
to those of figures 37 throu^ 4l. 

Figures 46 throu^ 48 illustrate the variation of lift in pure 
translation throu^out a cycle of motion. The amplitude of the sinusoidal 
translatory motion is O.9 inch in all cases. 

Figure 46 Illustrates the effect of increasing mean angle of attach 
for the airfoils in a pure translatory motion. It can be seen that the 
loops enclose very nearly equal areas except for the case -sdiere the air- 
foil is completely stalled. The change in size of these work loops for 

a given airfoil is dependent upon K = ^ in addition to the mean angle; 

hence the plots of work per cycle in translation (figs. 23 throu^ 25) 
reveal a more complete picture. The more general trend in translation 
is an Increase in negative work per cycle up to a mean angle approxi- 
mately equal to the static stalling angle and then a sharp decrease in 
negative work per cycle for angles above the static stall. 

Figure 47 Illustrates the effect of increasing frequency at a con- 
steint mean angle of 10°. As in the case of pitching motion, the enclosed 
area increases in size with frequency. This seems to be a general tend- 
ency for airfoils in translational motion. 

Figure 48 shows congparatlve data for the three airfoils used in 
the tests. 

Figures 49 throu^ 56 were chosen as typical examples of the oscil- 
lograph data analyzed and presented in this section, and where possible 
these plates are actually the basic data which are presented as the 
variation of lift or moment throu^ a cycle of oscillation in figures 37 
throu gh 4l and 46 through 48. These figures are presented to show in 
a qualitative way the variation of the reactions with time and the order 
of magnitude of the extraneous reactions and oscillations encountered. 

The additional points described below will be evident after an examina- 
tion of the figures by the reader. 
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Figure kS throu^ 51 are for a harmonic pitching motion. In the 
four pictures of figure one can see the pure sine -wave at 0°^ a 
sli^t distortion of the sine vave at 10°, and at l8° a vave -which is 
periodic hut hears -very sli^t similarities to a sine -wave. 

Figure 50 illustrates the differences in stal3Jjig characteidstlcs 
and in recovery from s-fcall for the -three airfoils sho-wn in figure tt-. 

The records of figure 51 show the effect of increasing frequency 
on the shaip-nosed airfoil of figure . 

Figures 52 "throu^ 5^ are for a harmonic translatory motion. Fig- 
ure 52 shows the effect of increasing mean angle on the intermediate 
airfoil. Figure 53 shows differences in airfoils at a mean angle of 10° 
and a frequency of approximately 8 cycles per , second, -while figure 5^ 
shows the effect of increasing frequency on -fche sharp airfoil. It -will 
he noted -that in general the differences in corresponding records for 
translatory motion are much less -than for pitching motion. 

Figure 55 shows the record oh-bained from the hlunt airfoil at an 
angle very near -the static s -tailing angle. This record is not periodic, 
apparently indicating a condition of partial stall all "throu^ the cycle 
of motion -wl-th the amount of s-tall not directly related to Induced angle 
hut prohahly more a function of tunnel flow inconsistency. Figure 56 is 
ano-ther record -taken very near the s-tatic s-talling angle. It can he seen 
that the wing is stalled during some cycles and unstalled during others. 
The records of figures 55 uud 56 are indications of the poor degree of 
reproducihility of data at angles near -the s-tall as mentioned earlier in 
-this section. This situation occurs more frequently in translatory ' 
motion than in pitching motion, and its effects are more pronounced in 
translation. It is apparent that -the solution to -this problem is a 
tunnel with an absolutely uniform flow. 


Correlation with Ibq)erimen-tal Da-ta from Other Sources 

Other experlmen-tal da-ta which have been gathered on airfoils oscil- 
lating at hi^ angles of at-tack may he dl-vlded into -two general types; 

( 1 ) Critical speeds oh-tained in a manner similar to that used in conven- 
tional flutter tests, except -that -the airfoil is placed at an initially 
hi^ angle of at-tack in -the air streami and (2) tests -which measure the 
oscillating aerodynamic moment in pure pitch in a manner similar to that 
used in -the M.I.T. tests, hut -which neglect o-ther reactions and motions. 
The tests under -the category of item (l) above -were carried out for -the 
most part several years eigo. Several investigators have contributed to 
-these da-ta, Studer and Kaufman -whose -work is summarized by Yictory in 
reference 3^ Rauscher in reference 9^ and BoUay and Brown in reference 10 
being -the most noteworthy. 
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An even larger number of e^cperlments has been made in the category 
of item (2) above. Victory in reference 3 describes tests in -which only 
the out-of-phase component -was measured, and Bratt and Scruton in refer- 
ence 5 plot -tdie path of moment through a cycle of oscillation much in 
the same manner as is done in this report. 

Bratt and Wi^t in reference 7 have been able to measure bo-fch the 
in-phase and out-of-phase coni)onents of moment, thus making possible 
the confutation of vector anflitude and phase angle of the fundamental. 

One further interesting reference is available -tdiich falls into 
neither of the categories listed a^jove. Reference 11 by Bratt, ¥i^t, 
and Chinneck is an investigation of aerodynamic danfing in pure pitch 
as affected by initial angle and angular aoflitude of oscillation. 

High-angle- of r-attack flutter tests .- Rauscher and Studer -were among 
the first ejcperimenters to show -the sharp drop in critical speed as angle 
of at-tack is increased to-ward the stall. Studer ’s investigation showed, 
and since then other investigators have shown, that, for lew values of 
freq.uency peirameter, the drop in flutter speed may be as much as 6o per- 
cent of -the classical flutter speed. However it has also been shown 
that as frequency parameter is increased the percentage drop in flutter 
speed decreases so that at very hi^ values of K the classical and 
stalling flutter speeds tend to approach a common -value. 

Studer (reference 12) -was able to obtain some interesting results 
by varying the mass and geometrical properties of his wing. A minlTmim 
stalling flutter speed was found to occur when the elastic and gra-vlty 
axes were coincident so -that there -was no inertia coupling. A critical 
speed below even this minimum -was obtained by restraining the -wing from 
translatory motion and allowing only freedom in pitch. 

Studer also varied the bending torsion frequency ratio, and found 
that the critical speed in s-talling flutter was quite insensitive to 
such variations. In all of the cases investigated in Europe the flutter 
motion \ia.8 observed to be predominantly a pitching oscillation even 
thou^ the same physical configuration had primarily a translatory motion 
in classical flutter. However, BoUay and Bro-wn at Harvard observed a 
stalling flutter of a solid Duralumin propeller section in principally 
the first bending mode. The critical speed -was 127.5 miles per hour as 
con^iared -with a classical flutter speed of 900 miles per hour so that 
it is almost certain that the rather -violent -vibrations -were actually 
stalling flutter. It is possible that the first bending mode of the 
blade was equl-valent to a translation plus an in-phase rotation about 
an axis in front of or behind the blade, but -there are insufficient da-ta 
really to fix the mode of flutter in reference 10. BoUay and Brown 
differentiated between the usual s-tall flutter and the excl-tation 
associated -with Karman voiiices. Their results showed that, since -the 
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Karm^ vortex frequency -was directly proportional to velocity, moderate 
vltratlons built up at airspeeds ■wb.ere vortex frequency coincided with 
a natural frequency of the wing. These vibrations however were much 
less violent than the oscillations described above and the reduced 
frequency for Karman vortex excitation was about three times that for 
stall flutter. These same experimenters made a number of tests on 
another wing, and their results are used as a basis for some of the 
confutations and text In the next section. 

Measurement of aerodynamic reactions .- A large part of the data In 
reference 3 which Includes the work, of Bratt and Scruton (reference 5) 
has been oreplotted against the same variables used In the presentation 
of M.I.T. data. 

It can be seen that the M.I.T. and the British data cannot be made 
to coDfare quantitatively since there Is no theoretical correction such 
as that described In reference 4, which can be applied to account for 
changed pltchlng-axls position. Further the airfoil shapes used In 
obtaining the British data are different from the airfoils described 
in this report, and the Reynolds numbers are widely different. Also 
because of the low airspeeds used in the British tests considerably 
hl^er values of the frequency parameter were obtainable, so that the 
two sets of teat data do not cover the same range. 

Figures 57 throu^ 60 are taken from reference 3^ It is possible 
to obtain rather good qualitative confarisons with M.I.T. test data If M 
Is defined as the angular interval between the mean angle and the steady- 
state stalling angle of the airfoil (positive when la greater than 

the steady-state stall), and this parameter is used as a basis of com- 
parison. By using data from the M.I.T. tests for the low K range and 
the British data for the higher K lunge, and the hipest angle- of- attack 
range, figure 6l was constructed using parts of figures 20 throu^ 22 and 
57 throu^ 60. This figure shows approximately the way in which work 
per cycle (or aerodynamic moment) varies with K and a. First for low 
angles of attack the curves of work per cycle against K remain negative 
and near the theoretical curves, gireidually approaching zero work as Zsa 
approached zero. Rear Ai. equal to zero a positive-work area appears 
at very low values of K. As Ax is further increased the maximum value 
of positive work is increased and the range of positive work is narrowed 
but continues to move to hi^er values of K. 

Reference 3 has summarized the results of several Investigations 
Into the effect of Reynolds number on the out-of-phase cong)onent of aero- 
dynamic moment. The curves and tabulated data in this summary show that 
for hl^ Reynolds numbers (above 2 x 10^) and moderately hi^ values of 
frequency parameter (above 0.5) the effect of angle of attack on the out- 
of-phase con5)onent of aerodynamic moment is small. Thus in this range a 
pure pitching oscillation with a conventional rotational-axis position 
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would iDe danced. Bratt and Wi^t in reference 7 devote a considerable 
portion of their work to an Investigation of the effects of aspect ratio, 
which Is beyond the scope of this report. However these Investigators 
were able to measure both conponents of the osclllatoiy aerodynamic 
moment, so that their data can be plotted In a mnner similar to the 
M.I.T. data. Figure 62 Is a plot of work per cycle against K. These 
data conpare gulte favorably with M.I.T. data If the basis of conparlson 
Is Zs(x rather than ot^. 

Figure 63 Is a plot of oscillatory moment coefficients and phase 
angles and It can be seen that these curves are gulte similar to some 
of the data plotted from the M.I.T. tests. Theoretical curves are also 
plotted In figure 63 and phase angles particularly can be seen to agree 
q.ulte well with theory for low angles of attack. 

All of the data taken In Great Britain measures the oscillating 
aerodynamic moment less the stlU-alr Inertia and dasplng forces. In 
reference 7 coefficient to account for stlll-alr dapping has been 
obtained experlPientally, but none of the figures described above has 
such a correction applied. However In most cases the effect of this 
dapping Is small enough that the general trends of the data are preserved. 

Effect of oscillation amplitude .- In all of the British references 
It has been observed that the reactions on an airfoil oscillating at 
high angles of attack are functions of the anplltude of the oscillation 
as well as Reynolds number, frequency, velocity, and Initial angle of 
attack. 

Especially In reference 7 'was this particular phase of the problem 
Investigated. This series of experiments definitely established the fact 
that there are anplltude effects at hl^ Initial angles of attack, but 
for any given set of conditions only two anplltudes of oscillation were 
studied so that no curves of aerodynamic reactions against oscillation 
amplitude can be plotted. 

The authors of reference 11 attack the problem of amplitude effects 
In a dlffeient manner from the other experlPienters . Their thesis Is 
tha.t for any Initial angle of attack there Is some magnitude of angular 
deflection idilch can cause partial stalling of the flow and thus cause 
negative damping, of the motion. The Piechanlcal apparatus used In this 
series of experiments was an airfoil restrained against translational 
motion and spring-mounted In pltclilng so that. If It were deflected and 
released, a pure pitching oscillation resulted. At various values of 
reduced frequency, tests were run In which the airfoil was deflected 
through larger and larger angles from a given Initial angle of attack 
until an undamped pitching oscillation was observed on release of the 
airfoil. Thus a series of critical angular amplitudes was obtained, 
below which oscillations were damped and above which the oscillations 
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tended to diverge. Figure 6k is a .plot of such data taken directly from 
reference 11 . These curves show a stable range at low and at very hi^ 
angles of attack and a range of instability in between. A similar effect 
(with a constant angular an^ilitude) may be observed by referring to fig- 
ures 20, 21, and 22, and noting that with increasing initial angle of 
attack the curves move into the positive-work range and then back into 
the negative-work range. 


STALL-FLUTTEE AMLYSIS 


The term "stall flutter" is usually associated with a periodic self- 
excited vibration of an "airfoil" with definite mass and elastic prop- 
erties in a flow •sdiich is separated during part or all of the motion. ■ 

The basic difference between "classical flutter" and "stall flutter" is 
in the character of the flow. As pointed out by Victory in reference 3^ 
stall flutter is not likely to occur on conventional aiiplane wings, but 
it does occur on propeller blades and tvirbine blades and some basis for 
predicting stall flutter or at least some knowledge of its characteris- 
tics would be very helpful. Of course in these applications three- 
dimensional and cascading effects further con^Hdcate the picture. 

The experimental, evidence of stall flutter has-been obtained for 
the most paid; under controlled laboratory conditions -where the flow is 
t-wo-dimenslonal and, the airfoil is rigid and supported on springs. Even 
for these tests, howevep, no adeq.ua-te method of prediction has been 
evolved. 

This problem of prediction by analytical means can probably best be 
approached by examining the addltioral complexity introduced by the 
stalling of -the airfoil. Basically the mass and stiffness properties 
of the -wing are unchanged but a new -variable, -the initial angle of 
at-tack Oj^, has been Introduced into -tiie determination of the aerodynamic 
derivatives. Perhaps a better form for this new -variable is something 
like Zipc. = cx^ - cCs-fceLll ^ figures k2 through 45. Regardless of 

-the form, however, the problem becomes extiamely coD 5 )lex because the 
aerodynamic deri-vatives can no longer be sinply superiB 5 )osed as in 
classical- flut-ter -theory; that is, -the aerodynamic results of a pitching 
motion- cannot be separated from -that of a simul-taneous translatory motion. 
There is -the further effect -that even in a pure motion the Introduction 
of the new angle- of-attack -variable automatically requires consideration 
of the large effects of such factors as Reynolds number and airfoil shape. 

Thus an exact mathematical formulation of -the problem is out of the 
question at present and -the at-tack must be made along enqplrical lines. 
Unfortunately only a small amount of engjirlcal data has been published 
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althou^ a number of general characteristics of the usual stall flutter 
can be deduced from these tests. As the angle of attack is raised 
through the stall region and the flow fs stalled during larger and 
larger parts of the cycle of vibration: 

(1) The flutter speed drops severely 

(2) The flutter frequency rises slowly toward the torsion frequency 

of the wing 

(3) The motion appears to become predominantly pure torsion althou^ 

the axis of rotation is not in general the elastic axis 

(!*•) Usually the flutter speed, reaches a minimum and then rises again 
as the flow becomes completely stalled 

A study of these general characteristics has suggested the following 
plan of attack which is based on proceeding from the con^jaratively firm 
ground df classical flutter into the realm of stalled flows. 

(1) Attempt to take initial angle of attack into account by modifying 
the existing classical theory by 'vdiatever means suggest themselves through 
e:jq)erience or intuition. Preferably these analyses should automatically 
reduce to the classical case at low initial angles of attack. 

. (2) Deteimlne which general kinds of modification produce the trends 

which have been observed in actual cases of stall flutter. 

t 

(3) Attempt to correlate these flndJuigs with. existing aerodynamic- 
derivative data in the hope of discovering which aerodynamic terms play 
a dominant part in stall flutter. 

(i<-) Try to predict the possible variations of these critical terms 
using airfoil data ordinarily available. 

This plan of course is aimed at the kind of stall flutter in which 
the angle of attack does not exceed the steady- state stall angle by more 

than 8 ° or 10 °. 


Den Hartog’s Sinpliflcation 

Den Hartog (reference 6 ) shows that one type of transmission line 
galloping is actually a stall flutter where the single of attack of the 
ice- coated cable is of the older of For an analysis of this motion 

it would of course be foolish to try to modify classical aerodynamic terms 
rather than make an attenpt at estimtlng the aerodynamic derivatives in 
the conpletely stalled flow. Since the frequency of the galloping is 
only about 1 cycle per second (K O.Ol), Den Hartog, as a first attempt, 
works only with an assumed static or steady-state lift curve. Considering 
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the changes in angle of attack to he Induced hy the vertical motion, he 
shows that the motion is negatively dangied on the negative- slope portion 
of the lift curve. Althou^ an excellent qualitative discussion, it is 
suggested that if quantitative answers to this problem are desired the 
rotation of the transmission line c ann ot he overlooked. It seems certain 
that, with the very low torsional stiffness of the long spein of cable, 
the combined effect of the mass unbalance caused by the ice formation 
and the aerodynamic moment about the cable axis will produce changes in 
angle of attack of the same order of magnitude as those induced by the 
vertical motion. 


Vlctoiy’s Use of Experimental Data 

Turning to the more usual type of stall flutter. Victory in refer- 
ence 3 follows the first three steps of the outlined plan of attack to 
show the dominant effect of the torsional dan 5 >ing parameter. This 
parameter is equivalent to 


1 1/1 
4 2\2 


a) - (1 - a2)F - j 


G 

K 


in the usual expression for aerodynamic moment in the vortex-sheet 
theory: 



It had long been suspected by the British that the reduction in flutter 
speed near stalling angles of attack, was caused by a reduction of aero- 
dynamic torsional dan^jlng. Thus Victory decided to try modifying the 
classical theory by varying only the parameter 

Fortunately experimental measurements of i^^d been made at the 

National Physical Laboratories for an airfoil which was very similar to 
one used by Studer (reference 12) for actual stall- flutter tests under 
almost identical conditions. In each case the symmetrical airfoil had. 
a reference axis at one-half chord and was tested throu^ about the same 
Eeynolds number and frequency-parameter ranges. Thus, by this happy 
circumstance, it was possible to substitute measured values of into 
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ttie modified classical theory and then check the results against com- 
parable stall- flutter tests. Inasmuch as the agreement between calcu- 
lated and experimental flutter speeds is very good, an exami n ation of 
the procedure is probably well worth -while. 

The series of tests by Studer for which calculations -were made had 
a common value of torsional frequency and six: different bending frequen- 
cies so that the range of frequency ratio 'was from 1.06 to 2.12. 

Flutter calculations -were nade for these conditions for values of 
ranging from 0 to 1.0 so that plots of flutter speed, frequency, and 
frequency parameter = 2K against could be constructed. These 

plots have been reproduced as figure 65. A curve of "the classical 
variation of sttperimposed on figure 65(b) 

to show the amount of -Tarlation of torsional danping required at 
any value of the frequency parameter to produce flutter. The intersec- 
tion of this curve -with each of -the other lines Indicates the value of 
corresponding to the classical- flutter solution for each frequency 

ratio. These values of "ttLen determine the flutter speeds and 

frequencies as marked by circles in figure 65(a). 

In order to obtain the semiempirical stall- flutter points. Victory 
plotted the six theoretical curves of against ^2 same 

sheet as the N.P.L. (j^^^c) curves, -where they cut most of the experi- 
mental curves at least once. At .these points of intersection the cor- 
rect experimental correspondence exists between ^3 and for the 

particular angle of at-tack at -which the experimental curve -vras plotted. 

The flutter speeds corresponding to these -values of were then 

ob-tained from figure 65(a). The agreement bet-ween these calculated and 
Studer’ s experimen-tal flutter speeds is quite good, especially for the 
larger frequency ratios -where the flutter speed does not -vary so shaiply 
-with mentioned by Victory the calculated flutter speeds have 

their sharp drop at a sll^tly lower angle of at-tack than the experi- 
mental speeds largely because the static s-talling angle of Studer’ s -wing 
is about 2.5° b1 gher -than the N.P.L. wing. If the speeds had been 
plotted against Ax = - ttg-taii rather than o^, this difference -would 

disappear. 

Thus it appears that the decrease, althou^ not necessarily the 
disappearance, of torsional danplng is an important, almost predominant, 
factor in the stall flutter of this particular -wing. 

• 

In an attenpt to extend these findings to another range of parameters, 
some theoretical points similar to those of figure 65 have been calculated 
for wing II tested by BoUay and Brown (reference 10) and used by Mendelson 
(reference l) in his analytical approach. The results of these calculations 
appear in table X along -with results ste m m i ng from Mendelson’ s analysis. 
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Mendelson’s Modification of Classical Theory 


A some'what different approach to the stall- flutter prohlem "was 
en^Jloyed by Mendelson in reference 1 . Following to some extent the 
plan of attack outlined earlier in this section, the classical theory 
was modified in the stall range to take into account a suspected aero- 
dynamic hysteresis or lag. As stated by Mendelson, "The assunptlon was 
made that the absolute magnitude of the oscillatory aerodynamic forces 
and moments are the same at the stall as at zero angle of attack but 
that the vector magnitudes of these forces and moments are changed, this 
change being caused by the lag of aerod yn a m ic damping and restoring 
forces behind the velocities and displacements at stall, thus giving 
rise to a hysteresis effect." Denoting this lag by the angle 0 , the 
classical expressions for lift and moment become: 

icjfc _ i rt 4 - + 


4 qb 




KG + iKFe"^^^^ e 


- It. 
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Substituting these expressions in the usual bending- tors ion flutter 
relations 


+ Cblioe^“* - L = 0 
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expanding the determinant of the coefficients of hQ and and 

separating into a real and an Imaginary equation, there result two 
quadratics in (cDo/cd)^ with coefficients involving K and These 

may he written: 

If the properties of a given wing are substituted in and 

equation can he solved for K for a given value of ^ hy finding, the 
intersections of Fj_ and plotted against K. 

This was done hy Mendelson for the wing used hy BoUay and Brown 
(reference 10 ). For zero aerodynamic lag, that is, for 0 = 0 °, the 
solution gives the classical flutter speed of 79 miles per hour. For 
increasing values of 0 the flutter speed drops almost linearly to 
about 30 miles per hour at a 0 of 55°- Encouraged hy this result 
Mendelson proceeded to relate this change in 0 to the steady-state 
lift curve. 

A closer examination of his answer, however, would have produced 
disquieting results. The drop in flutter speed is also acconpanied hy 
a decrease in flutter frequency "vdiich is contrary to the Bollay and Brown 
tests, to the usual experimental results, and to Victory’s predictions. 

It also appears that the assumed values of aerodynamic lag 0 are, 
especially for moment due to pitch, opposite in sign to the experimentally 
observed behavior in forced motion discussed in the next section. For 
the assumed change in vector hut not absolute magnitudes, the only way 
for Mendelson’ s 0 to follow the trends indicated experimentally is to 
take on negative values, that is, values less than 360°. 

This second observation indicates that perhaps Mendelson was looking 
for the equality between Fq and F2 for the wrong range of 0. To 
check this observation a calculation was made at a reduced frequency 
of 0.30 which corresponds to a '0 of 35° on Mendelson’ s curve. As can 
be seen in figure 66, another intersection exists for this K at 0 = 351° 
which yields a hi^er flutter frequency but lower flutter speed than for 
the classical case. The frequency is somewhat higher than mi^t be 
expected as shown in table X but the trend is correct. 

Since observed stall- flutter oscillations are usually described as 
predominantly pure pitching jnotlons a check tos made of anplitude ratio 
and phasing between the motions for the classical case and for 0 = 23°, 
35°^ and 351°* A 0 of 23° corresponds to a calculated flutter speed 
equal to the experlmenteLl speed measured at an initial angle of attack 
of 11. 2°. Values of 0 of 35° and 351° correspond to a calculated 
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flutter K very nearly the same as the experimental K for the same 
test. The results as shown in table X indicate that for all three values 
of ^ the motion is predominantly pitch althou^ about quite different 
axes. For 3 = 10.7° (0 = 35°) the effective axis of rotation is about 
half a semichord ahead of the elastic axis; for 3 = l'J2.k° {f> = 351°) 
the effective axis is about the same distance behind the elastic axis. 

It is unfortunate that BoUay and Brown did not record either phasing 
or effective axis of rotation. 

A basic difference be Ween Mendels on* s anfl Victory's approaches is 
the number of terms in the classical expressions which are modified in 
the stall range. Victory modifies only one critical term, the torsional 
dan5)ing. Mendelson applies his factor to a majority of the aerodynamic 
terms. In order to con5)are the results of the two procedures calcula- 
tions were made at a reduced frequency of 0.3 for the Bollay and Brown 
wing using Victory’s modification. In addition calculations were made 
at the seune reduced frequency applying Mendelson* s "corrected” factor <f) 
only to moment in pitch. The results, as shown in table X, show that 
for this case the additional modifications in^josed by Mendelson had a 
veiy small effect, althou^ none of the anstrers check very well with 
the experimental values. 

This latter discrepeincy may indicate that another aerodynamic term 
should be modified in an as yet untried manner althou^ it is possible 
that small errors in wing properties may be greatly nagnlfied in the 
calculations because the frequency ratio is quite near unity. To inves- 
tigate which aerodynamic terms plby an important part in the critical 
balance at flutter, vector plots have been constructed to illustrate 
some of the solutions. As shown in figure 67 the plots are referred to 
a unit a vector with separate vectors for lift and moment due to pitch 
about the elastic axis and for lift and moment due to translation. Lift 
and moment are plotted in such a way that the con^jonents of the moment 
vectors perpendicular to a are a measure of the work per cycle in 
pitch and the conponents of the lift measured perpendicular to h multi- 
plied by the ariplitude ratio h/bag are a measure of the work per cycle 
in translation. 

It appears that in the cases- corresponding to stall flutter the 
energy balance is governed by the torsional damping, the lift due to 
pitch, the phasing between the lift due to pitch and the translational 
motion, and the emplitude ratio. Thus the aerodynamic reactions resulting 
from the translation are very small but the work per cycle associated 
with translation combined with lift due to pitch is a very inportant 
part of the energy balance. On this basis a more conplete evaluation 
of the available data on lift due to pitch seems to be desirable in 
order to determine what other modification of the classical theory is 
appropriate. Variation of torsional dauping is inportant, perhaps pre- 
dominant, but is not the conplete answer to stall- flutter prediction. 
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Haley's Prediction of Aerodynamic Moment 

Apparently, then, it is possible to predict to some extent the 
trends of stall flutter for a particular configuration if experimental 
oscillatory data, especially for moment in pitch, are available. The 
data, however, must be for the correct airfoil shape, Eeynolds number 
range, and reduced, frequency range. 

These data, of course, are hardly ever available and are very 
expensive to obtain experimentally. Thus a very necessary step in an 
adequate prediction technique concerns the deduction of this information 
from steady- state data normally available. These available data should 
contain insofar as possible the effects of airfoil shape and Eeynolds 
number. 

Mendelson (reference l) attempted to correlate his phase lag 0 
with the changes in slope of the curve of steady-state lift against 
angle of attack. This type of correlation may be a fair first approxi- 
mation for relating a correction to the theory to the characteristics 
of the airfoil section but it makes no allowance for the variation in 
0 with reduced frequency which definitely does occur. 

Haley^ attempts in the case of moment in pitch to use the curve of 
steady- state moment against angle of attack to predict the actual varia- 
tion of moment in the stall range including the effects of reduced fre- 
quency and anplitude of oscillation. Actually all of his work was done 
for one airspeed but, as will be shown, the method can be generalized 
quite easily. Some intelligent guessing and experience is needed to get 
reasonable and useful results by this method but the end result could 
qviite possibly be a satisfactory set t)f curves of torsional danplng 
against reduced frequency for a useful range of angles of attack. 

Since the original paper is not generally available, a detailed and 
someidiat improved development is presented here. 

If at a given airspeed and instantaneous angle of attack the actual 
value of moment coefficient (or moment) differs from the steady- state 
value by an amount a possible assumption is that the rate of return 

of Cjjj to the steady- state value is proportional to This can be 

written sinply: 


dCjj _ ACJ^J 
dt Tj. 


where Tj. can be thou^t of as a characteristic time. Equations of 
this form describe many slnple physical systems whose response to a 

Unpublished M.I.T. Master's thesis by H. C. Johnson and S. M. 
Heiley, "Analysis of the Aerodynamic Eeactions on Airfoils Oscillated 
in the Stall Eange," June 19^9- C-126 
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step function is a sin^jle exponential governed 'by the characteristic 
tiine or time constant T^. 

As mig^t he expected, however, this equation is a little too sin 5 >le 
to describe the variation of moment. 

In congparing the predictions of this eq.uation with the known 
behavior of moment during oscillations at angles of attack below the 
stall, it was found that is different for different frequencies 

of oscillation and that the very sin 5 )le definition of ACj^ is not 
adeq.uate. An arrangement similar to the static and dynamic load-line 
concept in electronics was found to be necessary for handJLing oscilla- 
tions at various initial angles and frequencies. Thus Niyi is now 
redefined. (fig. 68) as 


ACj^ = Cy - (2) 

where is the value of given by the steady-state characteristic 

for the Initial angle of attack a^, X is the "dynamic" slope of moment 
against angle of attack, and a is the instantaneous angle of attack 
referred to and is usually written a = Using the new 

definition of together with the basic differential equation 

S (1) 

dt Tp 

the following equations are obtained: 

d(A%) _ da 

dt T^ dt 


dt 






a 


Assuming in each case a constant frequency oscillation a = a^e 


lout 


ACm = - 


XcLo^r 
1 + 


(oiT^ + ±)e^^ 


(3) 
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Cm = Cm + ^ (1 - 

^ 1 + (olTr)^ 


The classical expression for moment due to pitch is (reference 2): 

Cm “ °Mi + I ' 2(2 + 2(1 - a2)KQ + (l + a2)K^ + 


* a)o + 2(i - a2)KF - (| - a)K 


Equating the two expressions for Cm^ 


oJDj. = - 


2(^ + a^G + 2^ - a^KF “ 

^(2 + a2^r 


X = I 2(| + a^F - 2^ - a^^KGt^ + a^K^ 1 + (0113,)^ 


Thus it appears that the sin^jle assumed differential equation is adequate 
for unstalled flow if oTj, and X are determined for a given reduced 
frequency from the above relations. 

The fact that the product olPp is -the iiiportant parameter rather 
them Tj. alone suggests that a better formulation is possible. The 
product offj, can be considered as a frequency ratio which conp)ares the 
frequency of the motion to a frequency which is characteristic of the 
aerodynamic moment at the given airspeed. By suitable juggling it can 
also be written that 


^ V b Kr 
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Thus there is a characteristic quantity Kj. 'which includes both Tj. 
and its associated velocity in such a fashion that is a function 

only of reduced frequency K since equation ( 5 ) becomes 


+ a)F - 2(1 - a2) 

KG + 

(1 

2(1 .. a)s + 2(1 - a‘ 

-)ke - 

(1 - 


(7) 


The basic equation (l) becomes: 


dCM dCj4 




(8) 


where 



■which shows that the rate of change of Qm based on a time scale related 
to the motion is proportional to \Aiere the factor of proportionality 

depends only on the reduced frequency K. Equation (4) becomes: 



•which corresponds to an ellipticad contour on the plane of moment against 
angle of attack. The center of the ellipse is at line 

Joining the center ■with and _ is of slope X; and both 

the eccentricity of the ellipse and the slope X depend only on the 
reduced frequency K. 

Proceeding now to oscillations -which are partly or altogether in 
the s'tall region, certain observations can be made from the experlmen-tal 
results which can be a guide in extending the analytical technique. 
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The variation of the moment with angle of attack apparently can 
assume quite a variety of shapes as shown in figures 69 and 37 throu^ 

Ifl hut certain general characteristics can he observed. For oscillation 
into and out of the stall region: 

(1) The unstalled portion of the moment variation is essentially a 
part of the ellipse easily determined from the equations based on the 
dynamic slope X. 

(2) In the stalled portion, the moment at any angle of attack 
appears to tend toward the steady- state value of moment at the angle in 
a manner very similar to that below the stall. Usually these stalled 
steady- state values can be approximated by a strai^t line with a 
decidedly negative slope. 

(3) The actual stall occurs at an angle of attack above the steady- 
state stall angle but usually below the peak angle of attack of the 
oscillation. (See figs. k 2 throu^ ^ 5 -) The recovery begins approxi- 
mately at the steady- state stall angle. 

For oscillations ^ich are entirely above the stall such as the curve 
on the right side of figure 69^ the moment variation is approximately 
elliptical and appears to be related to its steady-state slope in a 
manner similar to the relation below the stall. 


Following these observations, it is possible to make some assunp- 
tions -vdiich enable an atten^Jt at prediction of the moment variation. 

(1) For the unstalled portion of the cycle use the appropriate part 

of the ellipse given by equation (9) 

(2) Assume that the same basic equation 


dC; 


M 


d(cut) 




(8) 


holds in the stalled portion of the cycle 


(3) Assume that = 


VT, 


has the same value at all points of 


the cycle 


( 4 ) Eeplace the steady-state variation of moment against angle of 
attack rtth a line of slope A, in the unstalled range and one or more 
strai^t lines in the stalled range 
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Hie assumed constancy of Kj. can be checked either by direct 
measurement of slope on the experimental records of from the value of 
work per cycle obtained by graphical integration. An accurate measure 
of the slope, -while subject to large uncertainties, gives a direct 
measure of K^. at any point in the cycle, -whereas -the value obtained 
from work gives a sort of average for a cycle. Numerous measurements 
of bo-th types have been nade and the results plotted in figure 70 against 
the theoretical variation confuted for the correct value of a. The 
con5)arison reveals "the encouraging indication that Kj. for stalled flow 
is certainly the same order of magnitude as for unstalled flow. It 
should be remembered ho-wever that these measurements all come from data 
obtained at one velocity so that -the results are some-what limited in 
scope. 

A number of moment variations -irere conqputed using these basic 
assun5)tions by a computer -who -was familiar with the general shapes 
ob-bained experimentally but did not have prior access to the particular 
experimental plots corresponding to his predictions. In some cases it 
-was assumed that the angle at -which stall occurred was the peak angle 
reached; in o-thers, stall angles were inteipolated from the experimental 
cuzves of figures k -2 throu^ To make reasonable transitions at stall 

and recovery from stall from -the ellipse associated -with one strai^t- 
line section of the moment curve to the ellipse associated -with the other, 
step-by-step integration was used. The results are plotted in figures 71 
throu^ 73 can be conpared -with the corresponding experimental curves 
included in figures 37 -throu^ 4 l. A tabuleur summary of the results 
appears in -table XI. 

The predictions are reasonable approximations of the actual varia- 
tions especially in some very inportant respects. The work per cycle 
and over-all amplitude are in good agretement with the experimental values 
for all -the cases tried -which Included a considerable variation in 
frequency and airfoil shape. The phase angle, however, i 6 in poor agree- 
ment unless a good estimate is made of the stalling angle. 

Thus it appears that a reasonably reliable plot of work per cycle 
or deinping coefficient against reduced frequency could be con- 

structed for -various initial angles of attack and amplitudes of oscil- 
lation provided the pertinent steady- state moment data are aveillable. 

At present, of course, -the accuracy of such predictions is known to be 
reasonable only in the general Eeynolds number and reduced- frequency 
ranges for -which -the above calculations were carried out. 

As stated earlier in this section, it may also be desirable to be 
able to predict the variation in the stall range of lift caused by 
pitching. Althou^ no such attempt has been made, a procedure similar 
to that used above may be possible. 
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If It is asjumedj just as for moment, that 

dt " " Tl 


and 


ACl 



the expression for at low initial angles of attack becomes: 


+ 




1 + 


(«l)‘ 


(1 - iairL)e 


The classical expression for lift due to pitch is 


iojt 


(reference 2): 



Con^arison of the two expressions produces the relations 

K + 25 + 2ft - a)KF 
ouTl ^ 

+ 2 F - 2(1 - a^KG 


Xl = - + 2 F - 2(1 - a)K^ ji + (oITl)^ 

•which are similar to eq.uations (5) and (6). Of course, the product olEj^ 
can be replaced by a reduced- frequency ratio of the foim K/Kl so that 
Kjj and Xj^ are functions only of the reduced frequency. 

Thus lift due to pitch for low initial angles of attack can be 
described in a m a nn er very similar to that used pre-vLously for moment. 
Whe-fcher a similar procedure can be used at and above the stall is not 
yet kno'wn. 
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CONCLUSIONS 


The aerodynamic reactions on a wing oscillating sinusoidally at or 
near the stall vary In a periodic but definitely nonsinusoldal manner. 

The form of the periodic variations, -which Is usually evaluated In terms 
of Its fundamental con^jonents, depends on the mean angle of attack, "the 
frequency parameter, the airfoil shape, the angpUtude of oscillation, 
the Reynolds number, and. In the case of a pitching motion, on the loca- 
tion of the rotational aids. This large number of lD 5 >ortant variables 
Increases the difficulty of predicting, or even measuring, the aerodynamic 
reactions far beyond -tiie kno-wn conplexl-fcy of -fche classical case at zero 
mean angle of attack. 

Most of the measurements to date have been confined to observations 
of the aetodynamlc moment resulting from a sinusoidal pitching oscilla- 
tion, althou^ the M.I.T. data presented in this report are a valuable 
source of Infonnatlon on the behavior of lift as -well as moment In both 
pitching and translational oscillations. Representative results from 
•fche various sources of data have been reproduced In graphical form in 
this report but because of the large number of significant parameters 
involved only a fe-w general conclusions can be dra-vm. 

It appears that the -variation of torsional damping -with mean angle 
of attack, -which can be rather precipitous in the general reduced- 
frequency range of 0.1 to 0.6 for Reynolds numbers of the order of 10^ 
and less, becomes q.ulte small for con 5 )aratlvely large values of reduced 
frequency and Reynolds number. This is particularly significant in the 
ll^t of -fche demonstrated importance in stall flutter of variations in 
torsional dangplng. Airfoil shape affects primarily -fche suddenness and 
type of flow breakaway under dynamic stalling conditions and, of course, 
together -with Reynolds number determines the range of angles of attack 
in -which stalling will occur. This last observation illustrates the 
Inadequacy for this type of work of referring angle of attack to the 
zero- lift condition. A more useful and meaningful reference is the angle 
of stall on the steady-state lift curve at -fche test Reynolds number for 
the airfoil in question. The effects on torsional dan 5 )ing of amplitude 
of oscillation and of rotational-axis position, -while fairly pronounced 
in some instances, are for the most part q.uite thorou^ly obscured and 
distorted by the variations of the other parameters. 

No Ins-tances of negative damping have been observed in the case of 
pure translation althou^ it appears that negative damping may possibly 
occur for values of the reduced frequency of about O.O 5 or less. Inasmuch 
as stall flutter Is ordinarily predominantly a torsional motion and the 
aerodynamic reactions associated -wi-fch the rotational coniponent are of a 
much greater magnitude -fchan those associated -with the treunslatlonal part 
of the motion, it is likely that further sBeasurements for pure translation 
•will produce no useful results. 
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No serious atten^t has apparently "been made to predict the aero- 
dynamic reactions under dynamic stalling conditions except for moment 
in pure pitch as has been proposed by Haley earlier in this report. 

This rather promising techniq.ue vhich requires only the appropriate 
steady-state lift and moment curves Includes to some extent in its pre- 
diction the effect of all of the important parameters in -rfiat appear to 
be their critical ranges. Extension of this technique to the prediction 
of lift in pitch appears possible. 


Considerable data on time-average values of lift and^^moment are 
available for the first time in this report and should prove useful in 
estimation of average or meain deflections of an airfoil such as a 
propeller blade for -vdxlch a stall- flutter prediction is being atten 5 )ted. 

The obseivation and prediction of stall flutter Itself have led to 
several interesting conclusions. A.s mentioned above most cases of stall 
flutter have exhibited a predominantly torsional motion with critical 
speeds much lower than predicted or observed at zero mean angle of attack. 
This difference in critical speeds usually reaches a maximum for mean 
angle not far from the steady-state stall angle and becomes less and 
less pronounced as the reduced frequency becomes large. Stall flutter 
should be differentiated from the forced excitation near a natural fre- 
quency ■vhich is often associated -with the shedding of the -well-known 
Karmim vortex trail. Stall flutter involves ty j^teraction of airfoil 
motion and aerodynamic reactions -vhereas the Karman vortices are observed 
behind even stationary objects. 

The prediction of stall flutter has been approached by attempting 
to modify adequately the usual classical -theory and has been moderately 
successful in a limited range. All of the results so far obtained Indi- 
'cate the major ln^ortance of changes in torsional dan 5 )ing and -when 
adequate experimental information on this parameter -vras introduced into 
the calculations the predicted stall- flutter speeds approximated the 
experimental values quite closely. Ho-wever it also seems evident that 
even thou^ stall flutter is predominantly torsional in character it 
cannot be reduced to a slngle-degree-of- freedom problem. The motion 
appears to be torsional only because rotation and translaticm of the 
reference axis are either almost exactly in phase or exactly out of 
phase. In the several cases examined closely the critical condition 
involved a balance between the damping in the translational motion and 
the pitching motion, both referred to the elastic axis. The lift arising 
from the pitching motion collaborates -wi-th -vdiat seems like quite a small 
amplitude translational motion to produce a large amount of negative 
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damping. Thus althou^ the stall- flutter motion appears to he quite 
elementary, that appearance merely hides the basic cong)lexity of the 
energy balance. 


Massachusetts Institute of Technology 
Cambridge, Mass., October I5, 19^9 
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TABLE I 

EXEEHIMEimL VALUES OF MAGHITODES ADD EHASE ANGLES FOR BLUNT WING 
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TABLE II 

EXPEREMENTAL VALUES OF MAGNITUDES AND PHASE ANGLES FOR 

INTERMEDIATE WING IN PURE PITCH AND PORE TRANSLATION 

[velocity 95 niplij elastic axis at 37-percent chord; 

semichord, 0.464 ft; pitch amplitude, 6.08°; translation 
an 5 >litude, 0.9 in^] 



Pure pitch 

Pure translation 

(deg) 

K 

^LP 

Ihs 

^MP 

4iP 

K 

^LT . 

l^T 

Cmt 

hm 

0 

0.094 

0.271 

178 

0.0680 

335 

0.096 

0.0388 

281 

0.0085 

80 


.168 

.240 

185 

.0670 

339 

.166 

.0568 

277 

.0133 

71 


.246 

.235 

188 

.0665 

327 

.272 

.0820 

232 

.0268 

62 


.339 

.229 

198 

.0670 

320 

.327 

.1000 

250 

.0227 

67 

6 

.055 

.229 

180 

.0545 

351 

.095 

.0945 

258 

.0087 

78 


.094 

.229 

182 

.0670 

340 

.171 

.0990 

246 

.0122 

68 


.171 

.215 

185 

.0700 

341 

.253 

.1240 

242 

.0199 

66 


.253 

.209 

183 

.0670 

322 

.327 

.1570 

233 

.0288 

70 


.339 

.215 

183 

.0775 

336 






10 

.053 

.189 

175 

.0292 

35 

.088 

.0380 

281 

.0075 

129 


.094 

.217 

184 

.0305 

32 

.166 

.0674 

268 

. 0 U 4 

97 


.164 

.217 

180 

.0462 

343 

.252 

.0890 

266 

.0180 

64 


.253 

.219 

201 

.0502 

333 

.319 

.1070 

260 

.0232 

63 


.318 

.212 

189 

.0593 

SI9 






12 

.055 

.211 

178 

.0704 

• 59 







.091 

.219 

199 

.0656 

45. 







.156 

.234 

206 

.0522 

32 







.247 

.237 

204 

.0375 

31 







.314 

.216 

194 

,0445 

13 





»• 

14 

.096 

.288 


.0467 

65 

.090 

.0497 

334 

.0099 

158 


.172 

.332 


.0452 

4 i 

.165 

.1014 

307 

.0170 

138 


.251 

.335 


.0398 

338 

.250 

.1264 

283 

.0182 

128 


.340 

.267 


.0343 

329 

.336 

.1405 

261 

.0223 

113 

18 

.095 

.304 


.0620 

111 

.091 


Nonpe 

riodic 



.179 

.368 


.0775 

87 

.170 

.0787 


.0225 

271 


.252 

.346 


.0633 

4 i 

.248 

.1230 


.0294 

251 


.332 

.310 


.0734 

342 

.319 

.1811 


.0457 

218 

22 

.055 

.150 

247 

.045 

180 







.091 

.196 

218 

.048 

162 







.168 

.244 

166 

.077 

208 







.259 

.278 

157 

.093 

200 







.335 

.323 

212 

.108 

157 




- . 
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TABLE HI 

EXEERIMEriTAL VALDES OF MAGNITUDES AND EHASE AHQTJBS FOB SHARP WING 
IN FDBE PITCH AND PDHE THANSLATION 

[velocity 95 elastic axis at 37-percent chordj semichord, 0,48^ ft; 

pitch an^jlltude, 6.08°; translation an5>lltude, 0.9 in7] 


Pure pitch 


Pure translation 


0.095 

.169 

.26k 

.337 



Cu> 

1kS2 

Cmp 


0.2ll9 

182 

0.054 

351 

.23^ 

188' 

.057 

340 

.23k 

187 

.061 

329 

.22k 

195 

.078 

328 

.228 

176 

.038 

22 

.219 

181 

.044 

10 

.210 

185 

.051 

351 

.213 

199 

.068 

334 

.219 

184 

.0476 

52 

.23k 

180 

.0287 

37 

.255 

214 

.0615 

l4 

.255 

206 

.0663 

0 

.278 

202 

.0655 

350 

,1814- 

178 

.0447 

79 

.235 

188 

.0708 

58 

.3014- 


.0723 

38 

.3014- 


.0686 

18 

.306 


.09I40 

0 

.296 


.150 

6 

.179 

182 

.0609 

90 

.242 

189 

.0738 

79 

.330 


.1025 

59 

.322 


11092 

30 

.292 


.0844 

30 

.282 


.0546 

97 

,390 


.0587 

83 

.374 


.0529 

52 

.325 


.0538 

32 

.243 


.0448 

156 

.354 


.0577 

112 

.336 


.0670 

80 

.309 


.0816 

52 

.105 

237 

.025 

190 

.126 

226 

.030 

199 

.180 

218 

.047 

183 

.234 

212 

.064 

201 

-314 

212 

.108 

183 


0.089 

.165 

.2k6 

.3^- 


.0895 

.175 

.260 

.332 


Clt ^T ^MT 0MP 


0.035 196 0.0077 81 ^ 

.056 261 .0119 69 

.083 259 .0167 52 

.118 271 .0219 59 

.070 270 .0061 112 

.097 2^ .0065 76 

.110 251 .0107 73 

.162 2lj-3 .0170 53 



Lodic 

.0638 2l4 

■ .095lt- 152 

.1190 129 

.0095 267 
286 



_KACA^ 
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I 

i 


EOSHDOmi. TIME-ATZSMIE TAUtES^ OT HAOKiroilSa IT»i HIS mci AKD njBl ZSA2laLmOC 
[jalooltjr, 9? (Jutlc ulii 37- O-farcest diard; Malcharaj O.ltd^F ft; pltA ai^altaila, 6 . 06 °; tnmlatlcii i^Utwle, 0.9 iaT| 



^&loa« vartf cr«s.-plottsd to olitaln th* data at ooaetut K'l iraosoteil In IlpixoM 3^ thraoih 36. 
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y 

WQHC CICIS in HJBE "PJSCE ABD HJBE TRAH3TAIEICIH 

[velocity, 9^ oq^; elaetlc axte at 37-Percent <^oid; eealchori., Itj translation. 

aspUtole, 0.9 pitch anpHtode, 6.03^ 



Pore pit^ 

Fore translation 

®1 

( 4 eg) 

Blpnt ving 

Intencedlate 

ving 

Sharp vLng 

Blunt v±Qs 

Intermediate 

vlng 

Sharp vlng 



ffet vorfc 


net voxk 


net iroric 


Net vorh 


Det ■voric 


Jlet voric 


E 

(ln.-Ib) 

E 

(IH.-11)) 

E 

(In. - 11 ) 

E 

(In.-lb) 


(!n.-lb) 


(tn.-lb) 

0 

0.093 

-1.33 

0.091 

-1.85 

0.093 

- 1.23 

0.092 

-3.88 

0.090 

- 1.18 

0.089 

-3.92 


.163 

-2.13 

.168 

-2.10 

.169 

-2.21 

.i£9 

-7.11 

.166 

-6.80 

.165 

-6.58 


.212 

- 2 .* 

.216 

-3.06 

.261 

-2.66 

.2I9 

-10.99 

.272 

- 9.35 

.216 

-9.33 


.331 

-3.48 

.339 

-3.98 

.337 

- 3.52 

.336 

-13.62 

.327 

- 11.10 

.311 

-12.33 

6 

.056 

-.71 

.055 

-.83 

.051 

-.35 




-8.91 


-8.90 


.093 

-1.62- 

.094 

-1.82 

.098 

-.75 

.091 

-8.26 

.095 

.092 


.170 

-2.12 

.171 

-2.81 

.171 

- 1.66 

.169 

-11.51 

.171 

-11.13 

.168 

-10.80 


.253 

- 2.93 

.253 

-3.27 

.231 

- 2 . 13 

.251 

-II.OI 

.253 

-15.00 

.251 

-12.50 


.311 

-3.20 

.339 

- 4.00 

.358 

-3.55 

.338 

- 17.20 

.327 

-16.17 

.316 

-16.50 

8 





.057 

0 












.091 

:33 












.168 

.99 












.2I9 

.17 








• 




.327 

-.66 







10 

.056 

1.42 

.033 

.08 

.056 

0 








.093 

-. 5 *^ 

.091 

-.86 

.093 

.12 

. .090 

- 1.09 

.088 

- 1.35 

.090 

- 2.15 


.172 

- 1.19 

.161 

-1.99 

.168 

1.22 

.167 

-6.37 

.166 

.7.91 

.175 

-7.93 


.251 

-1.75 

.253 

-2.72 

.251 

1.67 

• 253 

-10.80 

.252 

- 10.44 

.260 

-16.30 


. 3>»9 

-2.50 

.318 

-I.30 

.336 

1.33 

.335 

- 11.16 

.319 

-12.18 

.332 

-21.50 




.362 

1.08 







12 

.057 

1-95 

.055 

1.9I 

.058 

0 








.092 

2.76 

.091 

2.60 

.091 

1.06 








.168 

2.67 

.156 

2-35 

.166 

2.26 








.218 

.89 

.217 

.59 

. 2 l 0 

3.88 








• 323 

-.19 

.311 

-.41 

.363 

2.76 







1 ^ 

.077 

1.58 












.090 

2.58 

.096 

2.08 

.095 

-.40 

.089 

Nonperiodic 

.090 

-3.19 

.091 

nonperiodic 


.166 

5.00 

.172 

3.28 

.171 

1.20 

.170 

-12.62 

.169 

-9.35 

.168 

- 2.10 


.250 

5.36 

.251 

1.00 

,252 

2.66 

.260 

-21.50 

.250 

-14.50 

.256 

-el.oo 


.331 

1.55 

.3I0 

-.50 

.338 

3.87 

.3I3 

-23.75 

.336 

-19.00 

.333 

- 28.60 








.089 

Noz^perlodlc 












stall®! 

-9.00 

-14.51 













-16.13 





16 

.050 

2.15 












.159 

1.78 












.241 

5.81 












.325 

5.93 











18 

.099 

.59 

.C95 

1.42 

.095 

-1.63 

-090 

Nca^erlodlc 

.091 

Noi^rlodlc 

.090 

-.37 


.171 

3.68 

.179 

2.76 

.170 

.39 

.168 

- 2.88 

.170 

-3.60 

.168 

-2.52 


.257 

1.27 

.252 

5.08 

.252 

2.29 

.257 

-5.63 

.246 

-9.00 

.252 

- 7.00 


.338 

8.05 

.332 

9.18 

.319 

1.80 

.3IO 

-12.89 

.319 

-19.30 

.331 

- 12,20 

22 

.055 

-.92 

.055 

-.19^ 

.058 

-l.lB 








.090 

-1.I3 

.091 

-2.09 

.093 

-1.69 








.167 

-1.83 

.168 

- 2.00 

.170 

-2.36 








.218 

- 2.14 

.259 

-1.35 

.257 

-3.16 








.337 

-2.85 

.335 

-1.57 

.337 

- 3.83 















XAE2 n 
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TAELE VII 

RESULTS OF HARMONIC ANALYSIS FOR INTERMEDIATE 

IN PURE PITCH- AT VARIOUS AIRSPEEDS 

{ai, l 4 °j elastic axis at 37 -percent cliord.; semi chord, 0.484 ft; 
pitch an^plitude, 6. 08^ 


Velocity 

■tr 

% 

Mi 

% 


n 


(ng)h) 


(in.-lh) 

(in.-lh) 

(in.-lh) 

(in.-lh) 

OmP 

(deg) 

65 

0.079 

1.70 

5.0 

-5.6 

7.5 

0.0302 

138 


.133 

2.24 

6.7 

- 4.7 

8.2 

.0331 

125 


.244 

3.02 

9.1 

.9 

9.1 

.0367 

85 


.366 

2.09 

6.3 

6.8 

9.3 

.0375 

43 


.484 

-.35 

-.1 

9.5 

9.6 

.0387 

354 


.550 

-.82 

-.3 

7.1 

7.5 

.0302 

341 

80 

.064 

2.25 

6.8 

- 11.1 

13.0 

.0525 

149 


.107 

3.14 

9.4 

-9.4 

13.3 

.0536 

135 


.202 

4.26 

12.8 

-3.4 

13.2 

.0532 

105 


.298 

3.76 

11.3 

4.3 

12.5 

.0505 

69 


.391 

2.78 

8.4 

11.3 

l4.o 

.0565 

36 


.455 

2.62 

7.9 

13.3 

15.4 

.0621 

31 

95 

.054 

2.02 

6.0 

- 14.2 

15.5 

.0625 

157 


.090 

3.04 

9.1 

-9.4 

13.3 

.0536 

137 


.166 

5.01 

15.1 

-5.6 

16.1 

.0650 

no 


.250 

4.11 

12.3 

4.9 

13.2 

.0533 

68 


.327 

4.58 

13.7 

13.3 

19.1 

.0770 

46 


.380 

2.14 

6.4 

18.2 

19.2 

.0775 

19 

100 

.048 

CM 

CO 

• 

2.4 

-17.5 

17.7 

.0714 

172 


.079 

1.88 

5.6 

-17.5 

18.5 

.0745 

162 


.148 

4.68 

13.9 

-12.3 

18.6 

.0750 

131 


.216 

5.90 

17.5 

-3.4 

17.8 

.0717 

101 


.289 

■ 4.75 

l 4 .l 

9.3 

17.0 

.0685 

57 


.328 

3.61 

10.7 

18.2 

21.1 

.0850 

31 




TABLE VIII 


THEORETICAL VALUES OF AMPLITUDES^ PHASE ANGLES, AND RET WORK 


PER CYCLE FOE HIRE PITCH AND PURE TRANSLATION 


[pure 


tranalatlon h, O.9 la. J pure pitch, 


6,08^ 


Cmt 

^LT 

Pm: 

^T 

Net -vrork 

K 

Cmp 

OlP 



Ret work 

Cycle 
(in. -lb) 

Cycle 
(in. -lb) 

0 

0 

90 

270 

0 

0 

0.0805 

0.335 

3 ^ 

180 

0 

.00117 

.00486 

86.87, 

267.48 

-.584 

.01 

.0793 

.330 

356. 70 

177.90 

-.378 

. oo 44 i . 

.0182 

81.54, 

264 -. 10 

-2.18 

.04 

.0760 

.313 

350.78 

175.83 

-.995 

. 00639 

.0264 

78.90 

262.80 

-3.14 

.06 

.0743 

.302 

347.66 

175 . 42 

-1.180 

.00828 

.0339 

76 . 4-1 

262. ol; 

- 4.03 

.08 

.0728 

.293 

344.94 

175.50 

-1.581 

.'01009 

.o 4 o 9 

7^.78 

261. 6if 

- 4.86 

.10 

.0717 

.284 

342 . 52 

176.00 

-1.790 

.01189 

. 04-77 

73.08 

261.52 

-5.67 

.12 

.0710 

.2965 

340.30 

176. 76 

-1.990 

. 01511 

.0600 

70.22 

261.96 

-7.15 

.16 

.0700 

.2625 

338.37 

178.81 

-2.295 

. 01838 

. 0714 - 

67.80 

263.05 

-8.53 

.20 

.0699 

.2520 

332.93 

181.68 

-2.62 

.0216 

.0820 

65.65 

264-. 56 

-9.82 

.24 

.0705 

.2445 

329.87 

184.83 

-2.92 

.0262 

.0975 

62.80 

267.48 

-n.72 

.30 

.0722 

.2365 

325.88 

190.00 

-3.34 

.029k 

.1072 

61.03 

269.66 

-12.90 

.34 

.0738 

.2335 

323.56 

193.57 

-3.61 

.03k2 

.1220 

58.58 

273.21 

-14.70 

.40 

.0772 

.2320 

320.55 

198.97 

- 4.08 
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lABLE DC 


>5 


nKOKEHiaAL AMOIi OP AmCK AMD ABGDLAB THLOCIK OP A 3 CRPDIL AT STALL IM A FORE PilUUMU HOTION 


4 

Blunt vlngj staady-atats o^jti - 15 ° 

Inteimedlate ylng} stoady-etate 


Sliarp Tdngj steady-state ttotall “ 

<u 

(deg) 


^taU “1 

8a 

“Wu 

J 

£fl,rtall “1 

di 

'^atall 


dOwH-an a'DOTe 


■^tall 


(dag) 

( radians/ B«c) 


(deg) 

(radiana/sec) 

(deg) 


(radlans/sec) 

6 

0.056 

Ho Btall 



0.055 

Ho stall 



0.054 

5.56 

0.56 

0.62 


.093 

T *Tr» 

M ^ 



.094 

T *7n 

-do- 



. 0 ^ 
T Tl 

5.90 

ITa \ 

.90 

.54 




M ^ 





do 



.253 

.339 

.055 

'1 '' 



.254 

.398 

.056 

tl T1 




, 3 Xh 

,056 



n 



fi 



ax) 

5.15 

0.15 

0.87 

4.92 

0.92 

0.95 

i.i 4 

.14 

1.66 


.093 

5.25 

.25 

1.36 

.094 

5.12 

1.12 

1.48 

.093 

3.20 

2.20 

2.40 


.172 

Ho atall 



.164 

Mo stall 



.168 

4.75 

^.75 

3.09 






‘<35 

.318 




.iCJj. 

.336 

.362 

5.70 

Ho stall 
Ho stall 

trTO 

P. ^ 4 . 


.3^ 




3. D 













12 

.057 

5.70 

2.70 

.53 

.055 

3.06 

1.06 

1.42 

.058 

0 

1.00 

1.75 ■■ 


.092 

5.97 

2.97 

.28 

.091 

4.24 

2.24 


.091 

1.17 

2.17 

2.69 


.184 

5.99 

2.99 

.32 

.156 

5.06 

3.06 

2.^ 

.167 

3.14 

4.14 

4.28 


.247 

6.00 

3.00 

0 

.247 

5.16 

3.16 

3.78 

.240 

5.22 

6.22 

3.56 


.323 

Mo Btall 



.314 

5.85 

3.85 

2.10 

.362 

5.19 

6.19 

5.50 


.057 

.84 

-.16 

1.70 



2.54 







.0^ 

1.95 

.99 

2.56 

.096 

2.54 

2.62 

.0^ 

-.68 

2.32 

2.84 


.166 

4.53 

iS 

3.26 

.172 

4.52 

4.52 

3.40 

.174 

.78 

3 . 7 « 

5.20 


• 2V) 

■ 5.48 

S.06 

.251 

Ho stall 



.252 

3.68 

6.68 

6.00 


.331 

Mo stall 



.340 

Ho stall 



.338 

5.11 

8.11 

5.33 

18 

.099 

- 2.49 

.51 

2.72 

.095 

-2.15 

1.85 

2.68 

.095 

-5.49 

1.51 

1.16 


.171 

.71 

3.71 

5.12 

.179 

.31 

4.31 

5.39 

.170 

-2.16 

4.84 

4 . 8 o 


.257 

3.57 

6.57 

6.23 

.252 

4.53 

8.53 

^96 

.252 

1.02 

8.02 

7.53 


.338 

5.66 

8.66 

3.35 

.33s 

No stall 



.349 

3.72 

10.72 

8.25 
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TABLE X 


! ELimER COHDITIOITS FOR BOLLAI AMD BROWN, WHTG II 

i jwing properties: a = -0.29 ; Td = 3*375 in. J % ” 0.228; = 0.5 (asaumed); 

' = 16 i. 2 ; cdo, = 87.2 radians/sec; usb “ 90.3 radians/sec; 1%/% = 1.08^ 

L . ^ 

i 



(deg) 

/c 

(n 5 )h) 

K 

(rsidlana/sec) 

(ieg) 

ho/b 

^3 

Experimental: 

oi >= -1.36° 


78 

0.19 

77.2 




ai => 11.2° 


57 

.30 

82.8 . 




Calculated: 

Classical 

0 

78.5 

.20 

81.8 

125 

0.159 

0.371 

Mendels on 

23 

57 

.245 

72.8 

18.6 

.402 


Mendelson 

35 

45.6 

.30 

71.4 

10.7 

.579 


Mendelson ("corrected") 

351 

64 

.30 

100 

172 

.384 

— — — — — 

^ applied to moment In 

333 

63.7 

.30 

99.4 




pitch, only- 

victory’s modification 


67 

.30 

105 

— 


.027 
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TABLE XI 


'J5 


BEBULTS EROM HALEY’S EMPIRICAL TBEOEY 



Wing 

Sharp 

Sharp 

Blunt 

Blunt 

Intermediate 

tti, deg 

12 

12 

10 

10 

■ 14 

CD, cycles per second 

7.64 

16,66 

4.25 

7.90 

7.90 

Stall at a max 

Yes 

Ho 

Yes 

Yes 

No 

W, experimental, in.-lL 

4.52 

5.85 

-1.07 

-2.98 

6.56 

W, analytical, in. -It 

4.48 

6.44 

-0.7^ 

-0,79 

7.39 

ji, experimental, deg 

135.0 

60.8 

347.6. 

322.0 

106. 0 

analytical, deg 

41.2 

80,8 

352.7 

352.0 

125.8 

Amplitude, experimental, in. -It 

19.4 

19.4 

15.0 

20.6 

21.6 

Anrplltude, analytical, in. -It 

20.2 

19.4 

17.5 

17.1 

27.7 
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Line shaft and driving V-belt sheave 
Driven sheave and driving gears 
Crank wheels 
Vertical crossheads 
Connecting rods 
Position indicator 
Steel driving bands 




Tumbuckles 

Adjusting spring for jiroviding 
tension in bands 


Electric motor and gear train used 
for adjusting spring tension 

Upper supporting structure 

Wind-tunnel top and bottom 


Leading- and trailing-edge fairing 
strips for end plates 

Lower supporting structure 

Phase-setting protractor 



(a) Oscillator mechaniBin. 


Figure 2.- Diagrammatic drawing of test setup 
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(1^ Selectively daa^ied 
accelerometers 


(1^ Vertical giiidirig 
(1^ Guide bearings 

Initial-angle adjusting UnV 
@ Airfoil 



(ti) Model and supporting mechanism. 
Figure 2.- Concluded. 
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Figure 4 .- Airfoil section and ordinates 




Figure 5*- Arrangement of instrumentation. 


Hi TOVH 
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Angle of attack, deg 


Figure 8 .- Static airfoil cliaracterl sties. Intermediate airfoil. 


Moment aDout 37 -percent chord, in.-lb/sq ft 




Lift, Ib/sq ft 
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.5 0 5 10 15 20 


Angle of attack, deg 

Figure 9 -- Static airfoil cliaracteri sties. Sharp airfoil. > 


Moment about 37 -percent chord, in.-lb/sq ft 





Station 

Crank angle j 9 
(deg) 

Angle of pitch, a 
(a « 6.08 sin 0) 
(deg) 

1 “A a or A h 

1 

Vertical position, 
(h ■ 0.95 sin 9) 
(In.) 

0 

90 

6.08 


0.900 

1 

108 

5.79 


.856 

A 

Z 

126 

4.92 


.728 

3 

144 

3.58 


.530 

4 

162 

1.87 


.278 

5 

180 

0 


0 

6 

198 

-1.87 


-.278 

7 

216 

-3.58 


-.530 

8 

,234 

-4.92 


-.728 

9 

252 

-5.79 

* 

-.856 

10 

270 . 

-6.08 

1 +A a or A h 

-.900 


Figure 11.- Method of graphical integration of vork per cycle. 
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0 .1 .2 

K 


("b) Intermediate wing 



0 .1 .2 , .3 

K 


(c) Sharp wing. 
Figure 13 •- Concluded. 






.e in pure pitch 
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(c) Sharp wing. 
Flgiire l6.- Concluded. 

































foot span 
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JLu 

1 

X 

X 

14 
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16 

5 

-o- 

la 

7 

r 

22 

11 


Figure 22.- Work per cycle in pure pitch. Sharp wing. 
V s: 95 miles per hour. 
























.2 


.3 


K 

(c) Sharp wing. 
Figure 30.- Concluded. 











Figure 33.- Time-average value of moment in pure pitch. 
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0 2 4 6 8 10 12 14 16 18 20 22 

Initial angle of attack, ai, deg 

(b) Intermediate v±ag. 



0 2 4 6 8 10 12 14 16 18 20 22 

Initial angle of attack, o^, deg 

(c) Shaip ving. 

Figure 33-- Concluded. 








MCA TN 2533 


97 




0.2 4 6 8 10 12 14 16 18 20 


Initial angle of attack, Oj^, deg 
(b) Intermediate wing. 




(c) Sharp wing. 
Figure 3^-- Concluded. 
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-.7 

-.6 

-.5 
-.4 

a 

o 

-.3 
-.2 
-.1 
.0 

0 2 4 6 8 >10 12 14 16 18 20 22 

Initial angle of attack, deg 



Initial angle of attack, deg 


(c) Sharp ving. 



Figure 35-- Concluded. 
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Initial angle of attack, deg 



Initial angle of attack, Oj^, deg 


(c) Sharp wing. 
Figure 36 .- Concluded. 










u 




(d) Oi = i4°; Ax = 0°; 

<JD = 7-9 cycles per second. 


Continued. 


o 

uo 




cbordy liu-lb yw -foot wsmn 


I 


I 

1 

I 

[ 

; 

, W 



(e) = l8° I /^sa » l|-^j (f) « 22°; Ax = 8°; 

CD = 8.2 cycles per second. o) = 7*71 cycles per second. 


Figure 37.- Concluded. 
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(a) (D = 4.33 cycles per second. (Td) oi =. 7. Sit- cycles per secoixd. 


Figure 38.- Eaperlmantal hyBteresls loops from moment in pure pltcli. 
Intermediate vlng. = 6°; / 5 <x = -8°. 


S 

VJ! 


at chard 






abotrt diced, in, -lb per foot «p«n 


iaglw of ettadc, d«g 


Abtflo qC attaolc, 


(a) (u ni 4.18 cycles per second, ("b) (s = 7*64 cycles per second. 

Figure 39 •- Experimental hysterosls loops from, moaent In pure pitch. 
Sharp vlng. <xj[ => 12°; /to = 1°. 


An^ of BttAok^ dog 


of att«oAj dof 



c) 00 =■ 11.0 cycles per second, (d) o) = 16.7 cycles per second 

Figure 39*- Concluded. 





obord, Iru-lb p«r foot tfm 


30 


i 





VO 


(a) 00 - 4.27 cycles per second, (b) O) = 7.72 cycles per second. 

Figure 1 | 0 .- Experimental hysteresis loops from moment In pure pitch. 
Sharp vlng. = 10 °; Aa = - 1 °. 



(c) a> •= 11.5 cycldB per second, (d) <D = 15.1|-2 cycles per second. 

Figure 40 .- Concluded. 
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Htnonb i ti - r - fttv efacrdj ixu^Xb pii* foot *p*^ 














stall ” steady^state stall 
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("b) = l 4 °. 

Figure 45 .- Continued. 











(a) Oi = 0°j Aa = -l 4 °; (t) 0^ = 10°; Ax = -4°; 

to = 7*69 cycles per second. (jo = 7*63 cycles per second. 

Figure k 6 .~ E 3 q)erimental hysteresis loops from lift in prure translation. 

Intermediate ving. 
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(c) Ct^ = Zto = 0°; (d) «« l8°j 

00 = 7-59 cycles per second. od » 7*77 cycles per second. 




Figure 46.- Concluded. 
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■1.0 -.s -.4 -J 


.2 .4 .« 1,0 

-l.O -,8 -.4 -.2 

,2 .4 .6 .8 1,0 



TrsnfllAtlon aj^iUtuldj tn» 


Trauslrtloii aoplltisdej In. 


s 


9 

_ 


(a) tD = U.U cycles per second. (Id) £D «» 8.00 cycles per second. 

Figure 47 .- Experimental hysteresis loops from lift in pure translation. 

3hEirp wing, cn = 10“ j = -1°. 
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lift Lra r.iOT.ent Kositxo.i 1 

1 — 
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1 

1 
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1 ■ ■ 

1 
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r 
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— rV1 
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V' 
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\ ^ 1' * ' 

/ ^1 
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f ! 

PI 



(ift Draj Moxent Position 



(a) Blunt ving. 
a = 10° j 

■ CD = 7-9 cycles per second. 
Lift attenuation, 1-5; 
moment attenuation, 1-5; 
drag attenuation, 1-5 • 


("b) Sharp wing, 
a = 10°; 

00 = 7.7 cycles per second. 
Lift attenuation, 1-7; 
moment attenuation, 1-10; 
drag attenuation, 1-7. 



i . » • ' 

1'"l. 

'r-l,./. . . 

- 

1 

. ‘ , j 

1 ' /A 

V/ '> 

nP V . 1 

' ' 1 

(c; Intermediate vlng. 

■V a = 10°; 

• N CD = 7.5 cycles per second. 

'1 

■' ■ 

jAPl 

■ . ■ >■ /: 

JoLX U a w UcJjLliti UXOxl^ ^ 

moment attenuation, 1 - 7 ; 
-H drag attenuation, 1 - 1 . 5. 

1 

r : 

■ , 1 ' 

T ^ 

/ 

-J 



Figure 50.- Effect of airfoil shape. Pure pitch. 
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a = 10°; 

0 ) = k.2J cycles per second, 
lift attemxation, 1-5; 
moment attenuation, 1-7; 
drag attenxiation, 1-5. 



a = 10°; 

05 = 7.72 cycles per second, 
lift attenuation, 1-7; 
moment attenuation, 1-10; 
drag attenuation, 1-7. 



(c) 


a = 10°; 

0) = 11.5 cycles per second. 
Lift attenuation, 1-7; 
moment attenuation, 10-1.5; 
drag attenuation, 1-7. 


Dscillatlon. Shai^ ving. Pure pitch. 
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(a) Blunt -wing, 
a = 10°j 

CD = 7 - 67 cycles per second 
Lift attenuation, 1-3 J 
moment attenmtion, l-2j 
drag attenuation, 1-1. 
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(c) Intermediate wing, 
a = 10°; 

CD = 7-62 cycles per second. 
Lift attenuation, 1-3; 
moment attenuation, 1-2; 
drag attenuation, 1-1. 


Figure 53.- Effect of airfoil shape. 



Pure tranelation. 



■f 


Li t 


Draf; ^ Position 


Moment 

- - 



Figure 55 •- Record taken near static stalling angle shoving nonperiodic 
reactions. Pure translation. Blunt ving. a = lij-®j OQ = ll-.06 cycles 
per second. Lift attenuation, 1-3 j moment attenuation, 1-T> drag 
attenuation, 1-3. 
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Figure 56 .- Record taken near static stalling angle shoving an erratic 
transition from stalled to unatalled flov. Pure translation. Blimt 
vlng. a = l4°; m = 11. 9 cycles per second. Lift attenuation, 1-7; 
moment attenuation, 1-7; drag attenuation, 1-5. 
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Figure 5T*- Work per cycle in pure pitch. Reference 3^ figure 2. 
%teady-Btate stall " elastic ails at 50-percent chord; 

Reynolds mnnber = 1.42 x 105. 
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(a) Eeyaolda number = 2.83 x 10^ ; elaBtic axis at 33-percent chord. 

Figure 58.- Work per cycle in pxire pitch. Reference 3, figure 4. 

“steady- state stall “ 
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0 .2 .4 

I 

(c) Elastic axis at 33-P®rcent chord. 

. Plgure 59 •- Con 
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Figure 59.- Work per cycle In pure pitch. Iteference 3, figure 15 
“pteady-Btate stall = ^'5 elastic axis at 33-pereent ehor-d; 

Eeynoldfl nvmiber = 1.^ X 105. 






Work per cycle^ in, -lb 



Figurb 66 . ~ Work per cycle In pure pitch. Reference 3, figure 16, 
Steady-state stall “ elastic axis at 33-percent chord. 
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Figure 62.- Work per cycle in pure pitch. Eeference figure 1. 
Steady-state stall = elastic axLs at ,33 -percent chord; 

Reynolds number = 1.42 x 10^. 
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Figure 63 .- Moment amplitude and phase angle in pure pitch, from harmonic 
analysis in reference J. %tall “ 13 °; = 6°j elastic axis at 

33 -percent chord; Reynolds number = 1.42 x io5. 
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Figure 65.- Effect of torsional damping coefficient flutter for 

different frequency ratios Reference 3- 
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Figure 67.- Vector plots of aerodynamic derivatives at flutter. 


Moment coefficient 












Moment, in. -lb 
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(a) = 12° j CD = 16.66 cycles per second. 

Figure 71* - Hysteresis loops from eiaplrical theory. Sharp airfoil. 


Moment, in. -lb 
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(a) cti = 10°; ca = 7-9 cycles per second. 

Figure 72.- i^ysteresls loops from empirical theory. Blunt Eiirfoll. 
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("b) OjL = 10°j 0 ) ■= k.2^ cycles per second. 
Figure 72.- Concluded. 
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Figure 73.- Hysteresis loops from eD5>irical theory. Intermediate airfoil. 
Oi = l 4 °; CD = 7.9 cycles per second. 
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